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ABSTRACT 

We present a detailed analysis of the radial distribution of dust properties 
in the SINGS sample, performed on a set of UV, IR and HI surface brightness 
profiles, combined with published molecular gas profiles and metallicity gradients. 
The internal extinction, derived from the TIR-to-FUV luminosity ratio, decreases 
with radius, and is larger in Sb-Sbc galaxies. The TIR-to-FUV ratio correlates 
with the UV spectral slope /3, following a sequence shifted to redder UV colors 
with respect to that of starbursts. The star formation history (SFH) is identified 
as the main driver of this departure. Both Ltir/I^fuv and (3 correlate well with 
metallicity, especially in moderately face-on galaxies. The relation shifts to redder 
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colors with increased scatter in more edge-on objects. By applying physical dust 
models to our radial SEDs, we have derived radial profiles of the total dust 
mass surface density, the fraction of the total dust mass contributed by PAHs 
and the intensity of the radiation field heating the grains. The dust profiles are 
exponential, their radial scale-length being constant from Sb to Sd galaxies (only 
~ 10% larger than the stellar scale-length). Many SO/a-Sab galaxies have central 
depressions in their dust radial distributions. The PAH abundance increases 
with metallicity for 12 -|- log(0/H) < 9, and at larger metallicities the trend 
flattens and even reverses, with the SFH being a plausible underlying driver 
for this behavior. The dust-to-gas ratio is also well correlated with metallicity 
and therefore decreases with galactocentric radius. Although most of the total 
emitted IR power (especially in the outer regions of disks) is contributed by 
dust grains heated by diffuse starlight with a similar intensity as the local Milky 
Way radiation field, a small amount of the dust mass (~ 1%) is required to be 
exposed to very intense starlight in order to reproduce the observed fluxes at 
24 yum, accounting for ~ 10% of the total integrated IR power. 

Subject headings: dust, extinction — galaxies: ISM — infrared: galaxies — ul- 
traviolet: galaxies 

1. Introduction 

Understanding the spatial distribution of interstellar dust is of particular importance 
for two reasons: it affects our view of galaxies at different wavelengths, by absorbing UV and 
optical light and reemitting it in the infrared, and it also constitutes an important element 
in the chemical evolution of the interstellar medium (ISM). As for the first issue, correcting 
for dust extinction is usually the main source of uncertainty when deriving properties such 
as the star formation rate (SFR), age or metallicity (Calzctti ct al. 1994; Buat & Xu 1996; 
Calzetti 2001; Perez-Gonzalez et al. 2003). This limitation not only applies to integrated 
data, but also to surface brightness profiles and color gradients. For instance, the so-called 
inside-out scenario for the formation of galactic disks predicts that the timescale of gas 
infall and conversion into stars increases with radius, leading to radial variations of the star 
formation history (SFH), which are in turn observationally translated into color gradients 
(de Jong 1996; Bell & de Jong 2000; MacArthur et al. 2004; Taylor et al. 2005). Given that 
the dust content also changes with radius, it is key to properly quantify that variation with 
direct tracers of dust extinction; otherwise, derived parameters such as the radial growth 
rate of disks might be biased (see Munoz-Mateos et al. 2007 and references therein). 
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While being an inconvenience when studying star formation, dust itself is also a key 
ingredient in the chemical enrichment of the ISM. Metals resulting from stellar nucleosyn- 
thesis are returned to the ISM, where they condense to form dust grains, some of which 
arc later destroyed and incorporated into new generations of stars. Elements that later con- 
dense to form dust grains are injected into the ISM at different rates (see e.g. Dwck ct al. 
2009 and references therein). For instance, carbon now locked up into Polycyclic Aromatic 
Hydrocarbons (PAHs) might have been originally produced in AGB stars, whose typical 
lifetimes are a few Gyr. Other elements now constituting larger grains may have been syn- 
thesized in more massive stars, which die as supernovae in shorter timescales. Thus, the 
PAH abundance and the dust-to-gas ratio are expected to vary with the age of the stellar 
populations and correlate with the metal (oxygen) abundance of the gas (see e.g. Galliano 
et al. 2008 and references therein). Note, however, that dust formation results from a long 
chain of poorly-understood processes, of which stellar lifetimes and yields are only the first 
link. Grain growth and destruction in the ISM must be also considered. Indeed, it is thought 
that only 10% of interstellar dust is directly formed in stellar sources, with the remaining 
90% being later condensed in the ISM (see Draine 2009 for a recent review on the subject). 

Addressing these important issues requires a multi-wavelength, multi-object spatially 
resolved analysis, which is now possible for nearby galaxies thanks to the Spitzer Infrared 
Nearby Galaxies Survey (SINGS, Kennicutt et al. 2003). The SINGS project has made 
use of the Spitzer Space Telescope (Werner et al. 2004) to collect IR data, as well as an- 
cillary data from other facilities, for a sample of 75 nearby galaxies, representative of the 
galaxy population in the local universe. Together with UV images from the Galaxy Evo- 
lution Explorer (GALEX, Martin et al. 2005) and HI maps from The HI Nearby Galaxies 
Survey (THINGS, Walter et al. 2008), this data-set provides the scientific community with 
unprecedented spectral and spatial coverage of the most representative nearby galaxies. 

In this paper we exploit the unique data-set resulting from the GALEX-SINGS-THINGS 
collaboration by obtaining multi-wavelength radial profiles for the 57 galaxies in the SINGS 
sample which are detected and resolved at the longest wavelengths. Our subsample still in- 
cludes galaxies representative of different morphological types (ellipticals, lenticulars, spirals 
and irregulars). The profiles themselves, along with other observational parameters such as 
asymptotic magnitudes, concentration indexes and asymmetries, can be found in the accom- 
panying paper Munoz-Mateos et al. (2009, Paper I hereafter). The present paper focuses on 
the radial distribution of dust properties. Finally, by comparing the profiles of the spirals in 
the sample with models for the chemical and spectro-photometric evolution of spiral galaxies 
(Boissier & Prantzos 1999), we are studying the radial variation of the star formation history 
in these galaxies in a self-consistent frame, where radial changes in the gas-infall and chemi- 
cal enrichment are considered. The attenuation profiles turn out to be essential to properly 
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accommodate the model predictions with the observed profiles. Ignoring extinction not only 
biases the overall flux level in the UV and optical bands, but also the different radial scale- 
lengths at these wavelengths, which are a direct result of the inside-out assembly of disks. 
The results of this comparison with theoretical models will be presented in a forthcoming 
paper and, combined with the present study of the dust properties, should help to better 
understand how the different components of spiral galaxies get assembled together. 

This paper is organized as follows. In Section [2] we briefly outline the main properties of 
the sample, describe the different data used in this work, and explain how the radial profiles 
were obtained. Section [3] deals with the radial distribution of internal extinction, and how 
it relates with other properties such as the UV color, inclination or metallicity. In Section H] 
we compare our profiles with the dust models of Draine & Li (2007) (DL07 hereafter), and 
analyze the radial variation of the PAH abundance, the dust mass and luminosity surface 
densities, the intensity of the heating radiation field and the dust-to-gas ratio. Finally, we 
summarize our main conclusions in Section [51 In Appendix |A] we provide some empirical 
relationships between the dust properties derived from the model and the observed flux 
densities at the near- to far-infrared bands. In Appendix [B] we discuss the possible systematic 
factors that might affect our results in deriving the parameters of the dust models of DL07. 



2. The sample, data and procedure 

The SINGS sample (Kennicutt et al. 2003) consists of 75 nearby galaxies which span 
the range in morphological type, luminosity and FIR/optical luminosity observed in the 
local universe. The SINGS galaxies were also selected to cover a reasonably wide range in 
other additional properties, like nuclear activity, spiral and bar structure, inclination, surface 
brightness and environment. Note, however, that no significant luminous or ultra-luminous 
infrared galaxy (that is, with Ljr > 10^^ Lq) is included the sample. The median distance of 
the SINGS galaxies is 10 Mpc, with all objects being closer than 30 Mpc. Since the sample 
is neither flux- nor volume-limited, its statistical power as a whole is limited. Nevertheless, 
the wealth of panchromatic data available for these galaxies, together with their proximity, 
makes it possible to carry out detailed studies on the physics of star formation at kpc scales, 
including the interplay between star formation and the ISM. 

Eighteen SINGS galaxies were not suitable for our purposes and were excluded from the 
present study: 

1. Seven galaxies, DDO 154, DDO 165, Holmberg IX, M81 Dw a, M81 Dw b, NGC 0584 
and NGC 4552, are not (or marginally) detected in at least one of the MIPS bands. 
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2. Nine galaxies, Mrk 33, NGC 1266, NGC 1377, NGC 1482, NGC 2798, NGC 3265, 
NGC 3773, NGC 5195 and NGC 7552, are unresolved in the MIPS bands. Thus, 
radial profiles for these galaxies would reflect the shape of the PSF rather than real 
structures. In NGC 7552 some spatially extended emission can be discerned, but it is 
too faint compared to the extremely bright nucleus. 

3. The only evident IR source within the optical extent of NGC 1404 at 70 and 160 fim 
is an off-center object at the northeast of the frame. Since it might not be related to 
NGC 1404, we also excluded this galaxy. 

4. The bright center of NGC 3034 saturates the MIPS 24 fj,m detector, thus precluding 
reliable photometry. 

Therefore, the final subsample consists of 57 galaxies, whose main properties are sum- 
marized in Table [TJ 

2.1. GALEX data 

Nearly all SINGS galaxies have been observed in the FUV (Ag// = 151.6 nm) and NUV 
{Xeff = 226.7nm) by GALEX (Martin et al. 2005) (see Tabled). The observations are 
performed simultaneously at both bands thanks to a beam splitter, but some galaxies lack 
FUV images since the corresponding detector had to be turned off due to intense solar activity 
or overcurrent events. The delivered images have a final pixel-scale of 1.5". Although the 
size of the PSF depends slightly on the position on the detector and the brightness of the 
source, the typical FWHM is 5-6". This resolution is similar to that of the MIPS 24 /xm 
images, and corresponds to a physical scale of ~300 pc at 10 Mpc, the median distance of 
the SINGS galaxies. The flux calibration is based on white dwarf standard stars. For the 
pipeline version used here (the same as in Gil de Paz et al. 2007), the estimated zero-point 
uncertainty is 0.15 mag. 

2.2. Spitzer data 

The reader is referred to Paper I for a more detailed description of the Spitzer data used 
here. Mid- and far-infrared observations of the SINGS sample were carried out using the 
Spitzer Space Telescope (Werner et al. 2004). The Infrared Array Camera (IRAC, Fazio et 
al. 2004) was used to image the SINGS galaxies at 3.6, 4.5, 5.8 and 8.0 /im. We made use of 
the images provided in the SINGS Fourth Data Delivery, which are based on the Version 13 
Basic Calibrated Data produced by the Spitzer Science Center. The delivered images have 
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a pixel scale of 0.75", and the FWHM of the PSF at each band are 1.7", 1.7", 1.9" and 2.0" 
respectively. These correspond to spatial scales of 80-100 pc at the median distance of the 
sample. The estimated photometric uncertainty is ~ 2% (Reach et al. 2005). However, the 
photometry needs to be corrected for the diffuse scattering of incoming photons throughout 
the IRAC array^, and these corrections have an associated uncertainty of ~ 10%. 

The Multi-band Imaging Photometer (MIPS, Rieke et al. 2004) was used in scan- 
mapping mode to observe the SINGS galaxies at 24, 70 and 160/im. The final frames are 
delivered with pixel scales of 1.5", 4.5" and 9.0", respectively, thus being integer multiples 
of the pixel-scale of the IRAC frames while still properly mapping the MIPS PSF. The cor- 
responding FHWM are 5.7", 16" and 38" at each band, probing physical scales of 0.28, 0.78 
and 1.84 kpc at lOMpc. The estimated zero-point errors are 4%, 5% and 12% at 24, 70 and 
160/im, respectively (Engelbracht et al. 2007; Gordon et al. 2007; Stansberry et al. 2007). 

2.3. HI data 

The HI Nearby Galaxy Survey (THINGS, Walter et al. 2008) used the Very Large Array 
(VLA) to map HI 21-cm line emission from 34 nearby {D < 15 Mpc) galaxies, most of which 
were also targets of SINGS and the GALEX Nearby Galaxies Survey. The observations 
were done using the B-, C- and D-array configurations. For details of data reduction and 
processing, see Walter et al. (2008). Here, we use profiles derived from natural- weighted 
moment-zero (integrated intensity) maps. These have a typical angular resolution of 11" and 
sensitivity to surface densities as low as 4 x 10^^ cm~^ once convolved to our 38" working 
resolutions (see below). These radial profiles were kindly provided by A. Leroy and F. 
Walter. Because THINGS includes data from the compact B array configuration, the maps 
comfortably recover extended structure in our sources. The column densities in the THINGS 
maps are estimated to be correct to within ±10%. 

2.4. Surface brightness profiles 

The reader is referred to Paper I for a more in-depth description of the procedure 
followed to obtain the surface brightness radial profiles; here we will just briefly describe the 
most important steps relevant to this paper. In our study of the spatial distribution of the 
dust properties, we are limited by the resolution of the MIPS 160 yum images (FWHM of 38"). 



http : //ssc . spitzer . caltech. edu/irac/calib/extcal/ 
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We convolved the GALEX, IRAC, MIPS and HI images with different kernels (see Gordon 
et al. 2008) in order to match the shapes and resolution of their PSFs to the MIPS 160 fim 
PSF. Prior to convolving the images, foreground stars, background galaxies and artifacts 
were masked and interpolated over to avoid contamination when degrading the images. 

The radial profiles were obtained using the IRAF^ task ellipse, measuring the mean 
intensity along elliptical isophotes with fixed ellipticity and position angle, equal to those of 
the fiB = 25 mag arcsec"^ isophote from the RC3 catalog^ (de Vaucouleurs et al. 1991). For 
those objects for which these parameters were not included in the RC3 catalog, we used the 
major and minor axis diameters and position angles available in NED. The centers of these 
elliptical isophotes were set at the coordinates shown in Table [TJ The semi-major axis of 
these ellipses were successively incremented by 48" (a step larger than the PSF FWHM), to a 
final radius at least 1.5 times the R25 radius (depending on the extension of each particular 
galaxy). While using radially- varying ellipticities and position angles is useful in detailed 
studies of galactic structure at a specific wavelength, a panchromatic analysis requires using 
the same set of fixed elliptical isophotes in all bands to measure the different fluxes in the 
same regions of each galaxy, and for that matter a fixed PA and ellipticity was found to be 
as appropriate as any other set of values obtained from a given band. 

Uncertainties in the surface photometry include the error of the mean intensity within 
each isophote, computed assuming Poisson statistics, and the uncertainty in the sky level. 
The latter comes from high spatial frequency errors (Poisson noise, pixel-to-pixel variations) 
and low spatial frequency ones (flat-fielding errors) (Gil de Paz & Madore 2005). 

The radial profiles were corrected for Galactic extinction as in Dale et al. (2007), using 
the color excesses from the maps of Schlegel et al. (1998) and the extinction curve of Li & 
Draine (2001), assuming Ry = 3.1. The final profiles are shown in Table [2J 

Note that our radial profiles should be taken with caution in very inclined galaxies. 
First, although the radial step along the major axis is larger than the FWHM of the PSF, 
this is not the case along the minor axis in edge-on or close to edge-on galaxies. Moreover, 
the outer regions might be contaminated by emission from the central ones when performing 
the azimuthal average. Finally, the observed UV and IR radiation in these systems will likely 
probe different spatial regions within each galaxy, due to the large amount of dust along the 
line of sight. 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 

■^Except for NGC 5194, whose original values were highly affected by its companion galaxy, NGC 5195. 
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3. Radial distribution of dust attenuation 

3.1. Radial extinction profiles 

As part of our analysis of the radial distribution of dust properties, we will first determine 
the radial variation of internal extinction, which is also necessary to recover the intrinsic 
radial profiles from FUV to NIR. These kinds of studies have been addressed by several 
authors following different methodologies. Boissier et al. (2004. 2005, 2007) used the radial 
change in the total-infrared (TIR) to far- ultraviolet (FUV) ratio to derive radial extinction 
profiles for nearby galaxies. A similar procedure was followed by Popescu et al. (2005) on a 
pixel-to-pixel basis for M 101. Prescott et al. (2007) obtained extinction profiles in Ha for 
the SINGS galaxies by comparing the Ha and 24/im fluxes of individual star-forming regions. 
The number of distant galaxies seen through a spiral disk can also provide an independent 
estimation of the extinction (e.g. Holwerda et al. 2005 and references therein). 

Here we follow the first method to infer the radial distribution of dust attenuation for 
the SINGS galaxies. Several studies (Buat & Xu 1996; Meurer et al. 1999; Gordon et al. 
2000; Witt & Gordon 2000; Buat et al. 2005) have shown that the TIR-to-ultraviolet ratio 
is a robust tracer of the internal extinction in star forming galaxies, in the sense that it 
depends weakly on details such as the relative geometry of stars and dust, the shape of the 
extinction curve, or the star formation history (SFH). Regarding this latter issue, Buat et 
al. (2005) only found significant deviations in systems with very quiescent SFH at present- 
day (that is, havig decayed with exponential time-scales < 2Gyr). In these galaxies, the 
general starlight radiation field would become an even more important source of dust-heating 
than it already is in normal star-forming galaxies. This precludes computing the internal 
extinction in ellipticals, lenticulars and also the bulges of spirals with the same calibration. 
Using synthetic SEDs of galaxies with different SFHs and attenuations, Cortese et al. (2008) 
further investigated the dependence of the TIR-to-FUV ratio on the mean age of the stellar 
populations, and confirmed that quiescent systems exhibit larger TIR-to-FUV ratios than 
more actively star-forming ones with the same extinction. To account for the extra dust 
heating contributed by older stars, they provide a SFH-dependent calibration to estimate 
the UV attenuation. Note that, as pointed out by these authors, this recipe should not 
be blindly applied to ellipticals, since their FUV emission does not seem to be linked with 
recent star formation activity. We include these galaxies here just for completeness, but any 
further interpretation of their observed TIR-to-FUV in terms of attenuation should be done 
with caution. 

Using the low (48"-step) resolution profiles at 8, 24, 70 and 160 /^m, we built TIR 
(3-1100 /xm) profiles using the weighted sum proposed by DL07, after having subtracted 
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the stellar emission at 8 and 24 yum (see Section [4.21) . This estimator of the total- infrared 
luminosity constitutes a slight improvement over that of Dale & Helou (2002) using the 
MIPS bands alone, to the degree that including the 8 fim flux seems to reduce the scatter 
associated with differences in PAH abundances^. 

We then measured low resolution profiles on the convolved GALEX FUV and NUV 
images and estimated Apuv and Anuv from the TIR-to-FUV and TIR-to-NUV ratios, re- 
spectively. The attenuation at both wavelengths was computed using both the fits of Buat 
et al. (2005) and those of Cortese et al. (2008). Note that at small subgalactic scales, ra- 
diative transfer could limit the usefulness of the TIR-to-UV ratio as an extinction tracer, in 
the sense that UV photons emerging from a given region could heat dust in another region. 
Considering that our radial sampling is relatively coarse and that we average our data within 
elliptical annuli, this effect —if present— should be small (except maybe in extreme objects 
with very different UV and IR distributions, which are not common in the SINGS sample) 

The extinction profiles we derive are shown in Table [31 In Fig. [1] we show all the FUV 
extinction profiles for the galaxies in the sample, normalized to the optical radius R25. A 
general trend of Apuv decreasing with radius is clearly seen. The NUV profiles are not 
shown, since they exhibit similar behavior as the FUV ones. 

The extinction in the top panels has been computed with the calibration of Buat et al. 
(2005), whereas in the bottom ones we have used the SFH-dependent recipe of Cortese et 
al. (2008). While the equation relating Ltir/I^fuv and Afuv is unique in the case of Buat 
et al. (2005), Cortese et al. (2008) parameterized the coefficients of such a conversion as a 
function of observable colors that depend on the overall SFH. We followed the prescriptions 
given by these authors and determined the particular conversion between Lttr/ Lpuy and 
ApuY for each annulus depending on its observed (FUV — H) color. The latter was estimated 
assuming that {H — 3.6 /im) ~ —1.03 which, according to the stellar populations models of 
Bruzual & Chariot (2003), is the typical color exhibited by star-forming galaxies, and it is 
quite independent of their SFH. 

In SO / a-Sab galaxies and the bulges of later-type galaxies, the attenuation derived from 
the age-dependent calibration is ~ 0.5 mag lower than the one yielded by the age- independent 



''We have checked that the particular choice of this recipe is not critical, since both estimators yield 
almost equal TIR luminosities (on average, the TIR values derived from the DL07 calibration are 5% larger, 
with a rms of 5%). The DL07 recipe is a proxy for the actual TIR luminosity derived from the DL07 
models, designed to be valid over a wide range of starlight intensities and PAH abundances. In Section [4] we 
will derive Ltir from the model fitting, along with the remaining model parameters. On average, for the 
particular range of starlight intensities and PAH abundances in our sample, the TIR luminosities yielded by 
the estimator are only 6% larger than those from the model, with a rms of 3%. 
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recipe. The difference is negligible in the disk-dominated regions. Given that the SINGS 
galaxies span a wide range in star formation histories, both among them and within different 
regions of the same object, the results discussed in the rest of the paper refer to the extinction 
derived following the method of Cortese et al. (2008), except when mentioned otherwise. 

The overall level of extinction varies along the Hubble sequence, reaching a maximum 
in Sb-Sbc galaxies. On average, in these spirals the attenuation in the FUV ranges from 
~ 2.5 mag in the central regions to ~ 1.5 mag in the outer ones, although with large scatter. 
The extinction is ~ 1 mag lower in earlier spirals, and it goes below 0.5 mag in Sdm spirals and 
irregulars. Note again that the attenuation in ellipticals and lenticulars is highly uncertain, 
since even the age-dependent recipes of Cortese et al. (2008) might fail in these galaxies. 

In order to quantify the attenuation radial gradients, we performed a linear fit to the 
^Fuv('") profiles, without including the bulges of spirals for the reasons mentioned above. 
This exclusion was done visually: bulges produce a steep central rise of the 3.6 /xm luminosity 
above the exponential disk, and the FUV emission is significantly reduced in the central 
regions. This yields a sharp change in the (FUV— 3.6 /im) that, together with a visual 
inspection of the image, can be used to roughly delimit the bulge- and disk-dominated 
regions of the profiles (see also Paper I). The results are shown in Fig. [2], where the gradients 
are also expressed in terms of the R25 radius and the radial exponential scale-length of 
the 3.6/im profiles, defined so that /3.6/im oc e~^^"^■'^^^^. The scale-length aa.e^um was 
computed by fitting the profiles measured on the convolved 3.6 fim images (again, after 
excluding the bulges). Since at this wavelength the luminosity traces the stellar mass, this 
can be considered as the stellar mass scalelength. Most galaxies exhibit negative attenuation 
gradients, and the dispersion is larger in spirals of intermediate types. 



3.2. The IRX-/5 relation in normal disks 

3.2.1. Estimating the extinction from UV data alone 

In spite of the importance of correcting UV data for internal extinction, when FIR 
data are not available it is not possible to apply such corrections following the methods 
described in Section 13.11 In this regard, the slope of the UV spectrum [3 —or, equivalently, 
the (FUV— NUV) color— has been proposed as an indirect tracer of dust attenuation in 
starburst galaxies (Calzetti et al. 1994; Heckman et al. 1995; Meurer et al.l995, 1999). 
While the infrared excess (IRX) relative to the UV seems to be tightly correlated with f3 
in starbursts, later studies (Bell 2002; Buat et al. 2005; Seibert et al. 2005; Cortese et al. 
2006; Gil de Paz et al. 2007; Dale et al. 2007) have shown that the so-called IRX-/? relation 
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shows rather large scatter in normal star-forming galaxies. A similar behavior has also been 
observed for star- forming regions within galaxies (Calzetti et al. 2005). The correlation shows 
a wider spread, and is also systematically shifted to redder (FUV— NUV) colors. Using radial 
profiles from GALEX and IRAS, Boissier et al. (2007) confirmed that the IRX-/3 relation 
differs between starburst and normal star-forming galaxies. The relation is clearly shifted 
with respect to that for starburst galaxies, but the dispersion in their case seems to be 
reduced when using radial profiles instead of integrated photometry, possibly due to not 
mixing different stellar populations from the bulge and the disk. 

Here we extend the work carried out in that study by deriving a new radially resolved 
IRX-/9 diagram using data from Spitzer, which has better angular resolution than IRAS, and 
also probes colder dust by reaching a bit further into the FIR. An analysis of the IRX-/3 
plot for the SINGS galaxies, derived from integrated photometry, can be found in Dale et 
al. (2007). 

Throughout this paper we assume the definition of the UV spectral slope Pglx given by 
Kong et al. (2004): 

/?GLX = ({^/^v) - log(/A,Nuy) ^ 2.201(FUV - NUV) - 2 (1) 
log(AFuv) - log(ANuv) 

In Fig. [3^ we show the IRX-/5 plot for each of the low resolution radial profiles for the 
SINGS galaxies, classified according to their morphological type. Each data-point represents 
a given radial bin, and those belonging to spiral galaxies are also color-coded according to 
the (FUV— 3.6 /im) color of that bin (corrected for internal extinction), which can be used as 
a measure of the specific star formation rate (SFR per unit of stellar mass, sSFR hereafter) 
or, equivalently, the present to past-averaged star formation rate, usually referred to as the 
birthrate parameter b (Scalo 1986). This color scheme cannot be applied to ellipticals, since 
their FUV flux is not necessarily linked to star formation (Burstein et al. 1988; O'Connell 
1999; Boselli et al. 2005). As for the irregulars, the relative distribution of dust and stars 
is usually patchy (and not necessarily axisymmetric), so their Ltir/Ivfuv profiles must be 
considered with caution. As a reference, the right vertical axis shows the extinction in the 
FUV derived from Lxm/i^Fuv with the fit of Buat et al. (2005), although the color-coding 
is based on the (FUV— 3.6yum) color corrected for attenuation using the SFH-dependent 
calibration of Cortese et al. (2008). We also show the mean IRX-/5 for the starburst galaxies 
studied by Meurer et al. (1999), as given by Kong et al. (2004), and the fit provided by 
Boissier et al. (2007) for the radial profiles of normal spirals. 

Although differences in SFH can introduce systematic deviations in this plot (see next 
section), it is desirable to get at least a rough estimate of the internal extinction for normal 
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spirals when FIR data are not available, as has been traditionally done for starbursts. Since 
we are interested in obtaining a fit valid for normal star-forming spirals, we have excluded 
three galaxies with intense starburst activity (NGC 4536, NGC 4631 and NGC 5713). More- 
over, we have only considered regions with (FUV — NUV) < 0.9, thus avoiding the large 
scatter in the TIR-to-FUV ratio at redder UV colors. We applied a non-linear least-squares 
algorithm to the remaining data-points, giving: 

LtIr/Lp^V = 100-30+l-15(FUV-NUV) _ ^ 54 (2) 

The resulting fit is shown in Fig. [3^, and allows for the recovery of log(LTiR/IvFuv) for 
normal star-forming galaxies with a residual rms uncertainty of ±0.27 dex. Instead of im- 
posing a UV-color condition, we could have used the extinction-corrected (FUV — 3.6 /im)corr 
color to exclude the most quiescent systems. However, such a criterion cannot be applied by 
an observer lacking FIR data, which are necessary to correct that color for internal extinction. 
While our UV condition includes a few 'red' systems [in terms of their (FUV — 3.6 /im)corr 
color], they lie well within the main relation delineated by the rest of points. The resulting fit 
is not significantly sensitive to these few points, and this selection on the observed UV color 
is more justified from a purely empirical point of view. Beyond (FUV — NUV) = 0.9 mag, 
the observed dispersion is too large to reliably estimate Ap^y using this fit. 

Note that our fit differs from the one derived by Boissier et al. (2007) from GALEX and 
IRAS data, the latter having a somewhat flatter slope at the reddest UV colors, possibly due 
to data-points embedded within the bulges (no conditions on the UV color were imposed 
when computing that flt). Nevertheless, both flts are in very good agreement in the regions 
least affected by contamination from older stellar populations. By comparing the TIR-to- 
FUV proflles of the galaxies we have in common with the sample of Boissier et al. (2007), we 
have checked that the agreement between the TIR-to-FUV ratios is excellent. On average, 
our values are 0.04 dex higher, with a scatter of 0.13 dex. 

In brief, in order to determine the attenuation in the FUV one must flrst estimate 
the log(LTiR/i^Fuv) ratio. If FIR data are available, they should be used to directly com- 
pute log(LTiR using, for instance, the calibrations of Dale & Helou (2002) or DL07. Once 
log(LTiR/i^Fuv) is known, it can be translated into Apuv using any of the recipes available 
in the literature (e.g. Buat et al. 2005). In order to prevent the attenuation in early-type 
spirals from being overestimated, one can rely on age-dependent calibrations such as that of 
Cortese et al. (2008), which requires additional constraints such as optical and/or near-IR 
measurements. In the absence of FIR data, Apuy can be estimated from the (FUV— NUV) 
colour via Eq. [2], but such a relation should be only employed in a statistical sense for rela- 
tively large samples of galaxies. Given the large scatter, the use of Eq. |2] is discouraged for 
single objects. 
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3.2.2. Dependence on the star formation history 

In Fig. [3^, only a few data-points follow the starbursts relation, and they indeed cor- 
respond to galaxies with intense star formation activity (e.g. NGC 4536, NGC 4631). Most 
galaxies, however, lie to the right in the diagram: that is, for the same amount of attenua- 
tion, they have redder UV colors than starbursts. Differences in the star formation history 
are the most likely explanation for this broadening (Kong et al. 2004, Calzetti et al. 2005). 
Besides the effects of dust, regions with more quiescent SFHs will be intrinsically redder in 
the UV due to their more evolved stellar populations. Indeed, there is a clear trend with 
the (FUV — 3.6 /im)corr color, which, as noted above, is a proxy for the present to past- 
averaged star formation rate. Note that the 'present' SFR, as derived from the FUV, is 
not instantaneous, but represents an average over the last ~ 100 Myr. Data-points with 
(FUV — 3.6/im)corr < 3 mag follow a reasonably well-defined sequence, parallel to the star- 
burst one; but redder regions (usually embedded within bulges) depart towards the zone 
populated by ellipticals. 

The three dotted curves in Fig. [3^ show the IRX-/3 relation predicted by Kong et al. 
(2004) for galaxies with different values of the birthrate parameter h. The comparison with 
(FUV — 3.6 /im)corr should be done with caution, since the birthrate parameter of the models 
is the instantaneous one. The empirical relation for starbursts closely follows the model 
predictions for 6 = 5, which falls in the range where starbursts are commonly found (6 > 2-3; 
see Brinchmann et al. 2004 and references therein). The bulk of our data-points are consistent 
with lower values of h. These model predictions for Apuv, however, should be considered just 
as average approximations, given that their uncertainties range from ±0.3 mag for h > 0.3 
to ±1 mag for lower values, owing to differences in the particular details of the SFH and the 
dust content (Kong et al. 2004). 

We further explore this trend in Fig. HJ where we have plotted the (FUV — 3.6 /im)corr 
color as a function of the perpendicular (i.e. shortest) distance from each data-point to the 
starburst relation, c/gtarb (see Kong et al. 2004). Although with considerable dispersion, the 
trend is rather evident: regions with intrinsically redder (FUV — 3.6 /im)corr colors —having 
then lower current-to-past star formation activity— are clearly located further away from 
the relation for starbursts. Interestingly, selection effects can blur or even make this trend 
vanish. For instance, if we only consider regions bluer than (FUV — NUV) ~ 0.7, no clear 
correlation between SFH and (istarb can be inferred. Therefore, one needs to explore a wide 
range in SFHs in order to see this trend. This could imply that although SFH seems to be 
driving this departure from the starburst relation, other factors such as dust geometry or the 
extinction law might be also contributing in different ways. It should be noted as well that if 
we choose another IRX-/3 curve as a reference instead of the starburst one (i.e. MW-type dust 
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with different geometries, for instance), all perperincular distances will change accordingly. 

This could explain the different conclusions reached by authors studying integrated 
properties of galaxies. Kong et al. (2004) found a correlation between (igtarb and different 
indicators of the SFH, like the i5„(4000) break and the Ha equivalent width. Cortese et al. 
(2006) estimated the birthrate parameter from Ha and if-band luminosities for an optically 
selected sample of normal star-forming galaxies, and found a weak correlation between b 
and dstarb, with considerable scatter. The fact that their galaxies belong to nearby clusters 
might be in part responsible for the observed scatter. Interactions with the the intra-cluster 
medium and the cluster potential well may likely affect their SFHs by removing gas from the 
disks and quenching their SF activity, thus progressively turning them into anemic spirals 
(see e.g. Boselli & Gavazzi 2006 and references therein). In their analysis of the integrated 
properties of the SINGS galaxies. Dale et al. (2007) found that most of the scatter towards 
redder UV colors in the IRX-/3 diagram was due to ellipticals and early-type spirals. Panuzzo 
et al. (2007) studied a UV selected sample of galaxies, and did not find any systematic 
deviation from the starburst relation that depended on b, computed from NUV and H- 
band luminosities. However, their UV selected sample did not contain objects redder than 
(FUV — NUV) ~ 0.7, thus making it difficult to infer any correlation, as Fig. H] demonstrates. 

In order to further illustrate these issues, in Fig. [3]d we have highlighted some tracks for 
particular galaxies, connecting the data-points of annular regions at different galactocentric 
distances within each particular galaxy. The innermost point of each profile is marked for 
clarity, and the integrated colors of each galaxy are shown with open stars. 'Normal' spirals 
follow their own IRX-/? relation, which is quite similar in shape to the average one (modulo 
global offsets in the overall UV color and/or extinction). In early type spirals with well 
defined bulges, such as NGC 3031 (M 81), the track clearly deviates towards redder colors 
as we move closer to the center, due to the contribution of more evolved stars in the bulge. 
This does not happen in late-type galaxies, like NGC 2403, given their smaller bulge-to-disk 
ratios, especially at UV wavelengths. 

The SINGS sample also includes some peculiar objects that do not follow these smooth 
trends. This is the case of NGC 4826, which is a clear example of an anemic spiral (van 
den Bergh 1976). Although it has a global radial extent of ~ 13kpc at 3.6 /im, the bulk of 
the star formation activity —as traced by the FUV or 24 fim images, for instance— seems to 
be restricted to the central 5 kpc. Therefore, as we move away from the center towards the 
outer and more quiescent regions, the track followed by this galaxy in the IRX-/3 diagram 
heads towards the region populated by bulges. Something similar happens with NGC 4569, 
another anemic galaxy in the Virgo cluster, where ram pressure stripping seems to have 
quenched the star formation in the outer regions (Boselli et al. 2006). 
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Systems hosting starburst activity, like NGC 4536, lie close to the empirical relation 
given by Kong et al. (2004). In very edge-on systems like NGC 4631 (which also happens 
to host starburst activity), geometry might also play an important role. The observed TIR- 
to-FUV ratio is probably larger than the one we would measure if the galaxy was face-on 
(besides the fact that this galaxy is probably more dust-rich than the average spiral galaxy) . 
However, we do not find any significant trend between the position in the IRX-/3 diagram 
and inclination, apart from this extreme case. At fixed metallicity, however, inclination does 
seem to play a role (see Section 1X51) . 

Kong et al. (2004) suggest that galaxies in their sample departing too much from the 
starburst relation require very quiescent SFHs, with very short time-scales of star-formation, 
not typical of spirals. They invoke an additional mechanism that can contribute to broaden- 
ing the IRX-/? relation, consisting of an extra burst of star formation at some point during 
the galaxy's lifetime, superimposed on top of an otherwise smooth SFH. While this can 
certainly add more scatter, we note that in our radial analysis most regions lying at large 
distances from the starburst relation are either bulges or regions of anemic spirals with clear 
signs of star-formation quenching. Indeed, there are 18 galaxies in our sample having points 
with (FUV — 3.6 /im)corr > 3. Five of them are ellipticals or lenticulars, and 11 are early-type 
spirals (SO/a-Sb), including the already mentioned anemic ones. Therefore, at least for these 
very quiescent systems, a rapidly declining SFH might be a more reasonable possibility. 

From the above analysis we can conclude that the star-formation history is possibly 
driving the departure of star-forming regions from the locus of starburst galaxies in the 
IRX-/5 diagram. However, the trends with different tracers of SFH (both here and in other 
studies) seem to be quite noisy and not always evident —they can actually disappear if the 
range of explored SFHs is not wide enough. Several reasons might explain this. First of 
all, observational errors might blur the offset between regions if their SFHs are not different 
enough. Secondly, realistic SFHs cannot be parametrized with a single quantity like b; 
indeed, the intermediate bursts proposed by Kong et al. (2004) were shown to contribute to 
the observed scatter. Besides, a number of additional factors such as the relative geometry 
of dust and stars, the shape of the internal extinction law and the IMF can also have a great 
impact (Burgarella et al. 2005, Panuzzo et al. 2007). 



3.3. Influence of metallicity and inclination 

Attenuation and metallicity are known to be correlated in starburst galaxies (Calzetti 
et al. 1994; Heckman et al. 1998), which can be interpreted in terms of increasing extinction 
at larger dust-to-gas ratios (see the discussion and references given in Section 14.3.61) . Such a 
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trend has been also observed in normal star-forming galaxies, both with integrated data (see 
e.g. Cortese et al. 2006) and with radial profiles (Boissier et al. 2007). These studies find an 
offset with respect to the starbursts relation, with normal galaxies being less attenuated at 
a given metallicity. 

In Fig. [5^ we study the radial trend between metallicity and the TIR-to-FUV ratio. 
Metallicity gradients for 22 SINGS galaxies were taken from the compilation of Moustakas 
et al. (2009, in preparation). These authors compute the oxygen abundances using two 
different calibrations: Kobulnicky & Kewley (2004) and Pilyugin & Thuan (2005). The 
Kobulnicky & Kewley (2004) values tend to be smaller than the Pilyugin & Thuan (2005) 
by 0.6 dex, on average. However, to first order the radial abundance gradients derived for 
these 22 SINGS galaxies is independent on the adopted calibration. In this paper we adopt 
the oxygen abundances derived with the Kobulnicky & Kewley (2004) calibration, and refer 
the reader to Moustakas et al. (2009, in preparation) for more details on this subject. 

Although our photometric profiles extend up to ~ 1.5 times the optical size of each 
galaxy, the HII regions where oxygen abundances were measured do not usually reach be- 
yond the optical radius, and are also missing in the central regions of some galaxies. In 
order to avoid unsafe extrapolations of the metallicity values, we only plot those data-points 
lying within the radial range determined by the innermost and outermost abundance mea- 
surements in each galaxy. The empirical fits of Cortese et al. (2006) for cluster galaxies 
and of Boissier et al. (2007) for nearby spirals are shown for reference, as well as the one of 
Heckman et al. (1998) for starburst galaxies. 

We confirm that the most metal-rich regions also suffer from larger attenuation. Despite 
the scatter, we can see that at fixed metallicity the more inclined galaxies appear to have 
larger TIR-to-FUV ratios: from our perspective, UV photons must travel through larger 
amounts of dust before reaching us than if the galaxy was to be observed face-on. We 
performed a linear fit to those points with b/a> 0.5, resulting in: 

12 + log(0/H) = 8.70 + 0.405 log(LTiR/LFuv) (3) 

In Section 13.2.21 we stated that no significant trend with inclination can be detected in 
the IRX-/3 diagram, yet the TIR-to-FUV ratio is clearly larger for nearly edge-on galaxies for 
a given metallicity. This apparent riddle is solved in Fig. [5)d, where we plot the metallicity 
as a function of the (FUV — NUV) color for the same elliptical annuli as in panel (a). Both 
magnitudes are clearly correlated, and inclination also plays an important role here, since 
the observed UV color is redder at larger inclinations for a given metallicity. The effect of 
inclination depicted in panels (a) and (b) is such that more inclined galaxies appear redder 
in the UV, but also have larger TIR-to-FUV ratios. As a consequence, decreasing b/a will 
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not significantly displace the corresponding data-points out of the main trend seen in Fig. [3l 
but would move them along it. 

Since the TIR-to-FUV ratio seems to be linearly correlated with metallicity for face-on 
galaxies, the distribution of filled circles in Fig. [5)d is analogous to the IRX-/5 plot. The solid 
line in that panel is obtained by combining Eqs. [3] and [2], not from a direct fit to the data in 
the plot. 

4. Comparison with physical dust models 

4.1. Description of the models 

In order to derive the physical properties of the dust beyond mere dust attenuation, we 
need to assume a given model for the chemical composition and size distribution of the dust 
grains, and take into account their spatial distribution relative to the heating stars. Modeling 
the spatial distribution of cold and warm dust, relative to stars, has been the subject of much 
attention, especially in edge on galaxies (see e.g. Xilouris et al. 1999; Popescu et al. 2000). 
In this paper we will employ the dust models of DL07, which constitute an update of those 
developed by Weingartner & Draine (2001) and Li & Draine (2001). These models were 
successfully applied to the integrated photometry of the SINGS galaxies (Draine et al. 2007, 
D07 hereafter). The reader is referred to D07 and DL07 for a more detailed description than 
is provided for context in the following paragraphs. 

The DL07 models describe the interstellar dust as a mixture of carbonaceous grains 
and amorphous silicate grains, whose size distributions are chosen to mimic the observed 
extinction law in the Milky Way (MW). Although other size distributions can be chosen to 
reproduce the extinction curves in the Large and Small Magellanic Clouds, D07 showed that 
they did not provide significantly better fits than the MW ones for the SINGS galaxies. 

The physical properties of the smallest carbonaceous grains in the models are those of 
PAH particles. The PAH abundance is characterized by the PAH index, Qpau, defined as 
the fraction of the dust mass in the form of PAH grains with less than 10^ carbon atoms. 
Models with different PAH abundances can be generated by changing the size distribution 
of the dust grains while keeping constant the overall average extinction. The fraction of 
neutral and ionized PAHs is the one estimated by Li & Draine (2001) for the diffuse ISM in 
the MW. The resulting spectra have fixed band ratios in the IRAC bands, but those values 
are representative of normal spiral galaxies. 

One also has to make several assumptions regarding the starlight heating the dust. In 
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the DL07 models, the specific energy density of the starhght is set to be a scaled version of 
the local interstellar radiation field of the Milky Way estimated by Mathis et al. (1983). The 
intensity of the incident starlight is then characterized by a dimensionless scale factor U. 

Not all the dust within a galaxy is exposed to the same range of starlight intensities 
[/, and the emitted SED will result from the superposition of the SEDs associated with the 
diffuse interstellar medium, star-forming regions, etc. Dale ct al. (2001) proposed a power- 
law description of the dust heating, and used this model to describe the global infrared SEDs 
of the SINGS galaxies in Dale et al. (2007). 

The DL07 models follow a slightly different approach. The majority of the dust is 
supposed to be located in the diffuse ISM, being heated by a radiation field with a constant 
intensity f/min- A smaller fraction 7 of the dust is exposed to starlight with intensities 
ranging from f/min to f/max? following the power-law expression mentioned above. This term 
is intended to represent the dust enclosed in photo-dissociation regions (PDRs), where the 
radiation field is much more intense than in the diffuse component. Therefore, the amount 
of dust dMdnst exposed to radiation intensities between U and U + dU can be expressed as: 

= (1 - 7)Mdust5(C/ - C/min) + TMdust ^ f/-" (4) 

min max 

D07 find that the precise value of a does not have a great impact on the quality of the 
fits, given that the relative contribution of the diffuse and PDR components can be now 
parametrized with 7, and suggest fixing a = 2, which works well for a wide range of galaxy 
properties. They also find that no particular value of f/max seems to be favored, and hence 
adopt a fixed value of f/max = 10^. In addition, they suggest using 0.7 < f/min < 25 when 
submillimetcr data arc not available, since smaller values correspond to dust temperatures 
below ~ 15 K, which cannot be probed with MIPS photometry alone. This lower cutoff for 
U prevents the fitting technique from invoking larger amounts of cold dust heated by weak 
starlight with very small values of C/min, although this would underestimate the total dust 
mass if large amounts of cold dust are actually present. Note, however, that for those SINGS 
galaxies with available submillimetcr data, D07 concluded that omitting these data does not 
seem to introduce any evident systematic offset in the derived total dust masses, although 
the scatter in the resulting dust masses may increase up to ~ 50%. 

For a given PAH abundance, the emitted SED results from the hnear combination of 
the SEDs associated with the diffuse and PDR component. Therefore, the shape of the SED 
is controlled by three parameters: gpAH, f/min and 7, and the dust mass surface density can 
be derived when normalizing the SED. 



Besides the models of DL07, the SINGS galaxies have been used as a benchmark for 
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other different dust models. For instance, da Cunha et al. (2008) developed a simple model 
to interpret the mid- and far-IR SEDs of galaxies consistently with the emission in the UV, 
optical and near-IR. While their prescriptions for the dust emission are simpler and more 
empirical than those of DL07, by linking the absortion of starlight from the UV to the near- 
IR with dust emission at mid- and far-IR, they are able to extract the properties of both 
stars and dust in a consistent way. The dust properties they derived for the SINGS galaxies 
are in good agreement with those presented in D07. 



4.2. Fitting procedure 

Using our IRAC and MIPS radial profiles, we can therefore estimate gpAHj f^minj 7 
and the dust mass surface density SMdust ^ function of the galactocentric distance, by 
finding the best-fit model for the SED within each elliptical annulus. Following DL07, by 
taking linear interpolations of the original models we built a grid of models with different 
PAH abundances (ranging from ^pah = 0.4% to 4.6% in steps of 0.1%) and values of 7 
(from 0% to 30% in steps of 0.1%). We looked for the best-fit model by minimizing the 
reduced computed over the IRAC 5.8 //m and 8.0 /xm channels and the three MIPS 
bands, after correcting them for stellar emission (see below). We used the quadratic sum of 
the photometric uncertainties and the zero-point errors as weights, plus an additional ~ 10% 
error at each band to account for the limited accuracy of the model, following DOT. The 1-a 
uncertainty for each parameter was derived by projecting the overall distribution over 
the one-dimensional space of that parameter, and then looking for the values which satisfied 
x' = Xmin + l (Press etal. 1992). 

The DL07 models provide the multi-walelength flux densities emitted by the dust com- 
ponents. Since the IRAC 3.6 /xm and 4.5 /xm bands are almost entirely dominated by stellar 
emission (see e.g. Perez-Gonzalez et al. 2006), we did not take them exphcitly into ac- 
count when performing the fits; however, we used the 3.6 /im fluxes to estimate the stellar 
fluxes at the other bands. We used the stellar populations synthesis models of Bruzual & 
Chariot (2003) to compute the flux ratios {Fi,)x/{F,^)3,q. For a fixed age of 13Gyr, and 
assuming exponential SFHs with timescales r ranging from 0.01 to 20Gyr, we obtained 
{F^)x/{F^)3s =0.660, 0.453, 0.269 and 0.032 for A =4.5, 5.8, 8.0 and 24/im, respectively, 
assuming solar metallicity. These values are close to the ones derived from a 5000 K black- 
body SED, and depend very weakly on r, with variations of less than 5% over the range 
of extinction-corrected (FUV— 3.6 /^m) colors of our data. Larger deviations can be seen 
for super-solar metallicities, and also for very young populations (age < IGyr), due to the 
contribution of AGB stars. However, the prescription given above yields good fits for our 
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SEDs. 

To check the possible systematic effects introduced by the stellar emission correction, we 
fitted the models twice for all galaxies, using our flux ratios and the ones estimated by Helou 
et al. (2004) using stellar population modeling from Starburst99 (Leitherer et al. 1999). 
Their scaling factors are somewhat bluer than ours; hence the derived PAH abundances 
are systematically larger, since the estimated stellar emission at the IRAC bands is smaller. 
Nevertheless, the median offset in the derived abundances is only ~ 5-10%. There is, however, 
a much larger difference in those regions where the PAH contribution to the 8 /zm flux density 
is less than half of the stellar one at that band (e.g. in bulges). The values of gpAH might then 
be unreliable in these cases, so we have excluded them from our analysis. As for the fraction 
7 of the dust mass exposed to high- intensity radiation fields, for 7 < 1% the resulting values 
using Starburst99 are ~ 10% smaller than the ones obtained when using the factors from the 
models of Bruzual & Chariot (2003), with no difference at larger values. Finally, the dust 
mass and luminosity surface densities, as well as the radiation field scale-factor U, do not 
vary noticeably. This is not surprising, since these parameters are mostly determined by the 
MIPS fluxes, where the stellar contamination is negligible. For the remainder of the analysis 
described here we adopt the factors obtained from the models of Bruzual & Chariot (2003). 

Further systematic effects might result from using a small number of bands when fitting 
the models. In Appendix [B] we discuss this issue in detail, and demonstrate that our results 
are consistent with the global values obtained by D07 using a more complete set of IR data. 

4.3. Results from the model 

4-3.1. Individual profiles and SEDs 

In Fig. E] we show the infrared SEDs over a range of galactocentric distances, along with 
the best-fitting models, for one of the sample galaxies: NGC 7331 (see the on-line edition 
for a complete figure-set for the whole sample). The corresponding numerical values of the 
model parameters are quoted in Table HI which is shown in its entirety in the online edition. 
Note that the results for NGC 7331 can be readily compared with those from Thilker et 
al. (2007). As explained above, each SED can be decomposed in two terms: a 'diffuse' 
component (dashed line), heated by a constant radiation field with U = t/min, and a PDR 
one (dash-dotted line), with < U < f/max- The latter is particularly important to 

account for the observed flux density at 24 /xm. The total resulting SED is shown with a 
solid line, and includes not only the emission contributed by dust, but also the one coming 
from stars (dotted line), estimated from the 3.6 /xm flux. The observed flux densities (open 
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rectangles) are shown along with the ones resulting when convolving the model with each 
band (filled circles). Although the 3.6 and 4.5 yum fiux densities are shown, they were not 
included when performing the fits, as explained in Section 14. 2[ 

The radial trends of the different model parameters are shown in Fig. [7] for each individ- 
ual galaxy (only NGC 7331 is included in the printed version). In all panels, the best-fitting 
values are shown with different lines, and the gray bands correspond to the estimated l-cr 
uncertainties. Panel (a) shows the radial variation of the PAH mass fraction, Qpah- Regions 
where gpAH is not reliable are identified with a dotted line instead of a solid one. In these 
regions, stellar emission at 8 fim is more than twice the one contributed by PAHs, so even 
though their fitting errors might be small, they are affected by large uncertainties resulting 
from the stellar subtraction. Therefore, these values will not be considered in subsequent 
sections. The physical scale of the MIPS FWHM at 160 /xm (38") at the distance of each 
galaxy is marked with a small horizontal line. 

In panel (b) we show the fraction of dust mass exposed to high intensity radiation fields, 
7. The dotted and dashed lines correspond to the fraction of the dust-luminosity contributed 
by dust grains in those regions, heated by starlight with U > U^am and U > 10^, respectively. 
Note the different scales for the mass and luminosity fractions: about 1% of the dust mass 
can be responsible for ~ 10% of the total dust luminosity. 

Panel (c) displays the radial trends of the minimum starlight intensity f/mim as well as 
an average value {U), which depends on both f/min and 7 (see Eq. 33 in D07). With this 
definition, we have Ldust oc (f/)Mdust- 

Finally, in panel (d) we show the radial profiles of both dust mass and luminosity surface 
densities, corrected for inclination using the semiaxis ratio in Table [TJ The open circles show 
the dust luminosity computed from the monochromatic fiuxes from 8 to 160 /im, using Eq. 22 
from DL07. These values compare well with the ones obtained by integrating the SED from 
the model (dashed line). 

4-3.2. Properties of the whole sample 

In Fig. [H]we show the radial variation of each parameter in the model for all galaxies in 
the sample, divided into bins of morphological type. The radial coordinate is normalized to 
the optical radius, R25. The median profiles in each panel are shown as a thick dashed line. 

The top row shows how the PAH abundance changes with the galactocentric distance. 
Note that we only plot those values of gpAH satisfying the fiux criteria defined in Section 14. 2| 
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nevertheless, the PAH abundances for the earhest types should be still considered with 
caution. The median profiles of Qpak in SO/a-Sbc galaxies look somewhat fiat, although 
some individual profiles exhibit a central depletion of PAHs, relative to the abundances at 
larger radii. The overall PAH abundance in these spirals ranges from 2% to 4.6% of the total 
dust mass. As we move towards later types, steeper radial gradients become more frequent, 
and gpAH decreases down to ~ 1.5% or even less for the latest types. As we will see later, this 
can be interpreted in terms of a connection between the PAH abundance and metallicity. 

The fraction of the dust mass that is heated by very intense starlight clearly decreases 
with radius. This component is mostly determined by the flux density at 24 fim, which 
usually outlines the spiral pattern in which star-forming regions are typically arranged. The 
spatial filling factor of the HH regions is thus larger at the innermost regions, where they 
swirl and crowd together, and where sometimes they coincide with rings, whereas at larger 
galactocentric distances they are more dispersed, so their contribution to the 24 /im surface 
brightness is not so important in comparison with the diffuse one. Indeed, the radial decrease 
in 7 is steeper in early-type grand-design spirals, while in late-type flocculent ones the 
gradient is flatter, since the spatial distribution of star-forming regions is more uniform. 
This is in agreement with the findings of Bendo et al. (2007), who also concluded that the 
24 /xm emission is more centrally concentrated in early- type spirals than in late-type ones. 
In any case, although the values of 7 are quite small (typically around 1%), such a small 
fraction of the total dust mass can be responsible for a non-neghgible fraction of the total 
emitted infrared power. 

The scale-factor Umm describing the intensity of the diffuse heating starlight also de- 
creases with radius. This quantity is affected by both young and old stellar populations, 
and the bulge can certainly be responsible for part of that heating, raising the central val- 
ues of C/min, especially in Sa-Sbc spirals. Additionally, the higher central concentration of 
star-forming regions in these galaxies may also increase the diffuse interstellar radiation field 
(besides boosting 7). In the outer regions of spiral disks, beyond r ~ 0.5 R25, the median 
value of ?7min is close to 1, meaning that the radiation field has the same intensity as the local 
MW one. Noteworthy, in the latest types both Ujam and 7 present higher median values, 
although with more dispersion. 

The radial variation of the dust mass surface density also depends on the Hubble type. 
The largest amounts of dust are found in Sb-Sd galaxies, where the surface density decreases 
exponentially from lO^MQkpc"^ in the central regions to lO^MQkpc"^ in the outer ones, 
at the R25 radius. Conversely, dust is less abundant in Sdm and irregulars, as well as in 
ellipticals and lenticulars, with SO/a-Sab spirals lying in between. 

It is worth noting that the profiles of dust mass surface density appear to be exponential. 
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and the radial scale-length seems to hold rather constant from Sab to Sd spirals, when 
expressed in units of the optical size. To quantify this behavior, we have applied linear 
fits to these profiles. We have excluded prominent bulges from the fit, in the same way as 
we did in Section 13.11 with the attenuation profiles. This is particularly important in many 
SO/a-Sab galaxies, in which the central regions usually deviate from the exponential trend, 
due to holes or depletions in the amount of dust (see the full version of Fig. [7] in the on-line 
edition of the journal). 

The resulting exponential scale- lengths adust, defined so that ^M^^^t oc e~^/"<*"=*, are 
shown in Fig. [9l When expressed in absolute physical units (panel a), the dust mass scale- 
length varies widely across the Hubble sequence from 1 to 7kpc. However, the distribution 
appears strongly peaked when normalizing by the optical size (b) and by the scale-length of 
the stellar mass profiles (c), the latter obtained from the profiles measured on the degraded 
IRAC 3.6 /im images, for consistency. The dust mass radial scale-length is equal to 0.2- 
0.3 R25 for Sab-Sd galaxies. The profiles are flatter in later types (although the dispersion 
is larger, as these systems might deviate from the disk-like geometry). In general, the 
dust scale-length seems to be ~ 10% larger than the stellar one. By applying a three- 
dimensional radiative transfer model to seven edge-on spirals, Xilouris et al. (1999) concluded 
that the dust radial scale-length is ~ 40% larger than the stellar one. The presence of a 
cold dust component (T < 15 K) with a larger scale-length could explain this difference, as 
this component would be detected in absorption by the radiative transfer modelling, but it 
would be unconstrained by our MIPS profiles. Conversely, some of the necessary geometrical 
simplifications on which this kind of 3D models rely could introduce some bias. For instance, 
the dust profiles are assumed to be exponential throughout the whole disk, which does not 
seem to be the case in many early-type spirals exhibiting central plateaus or holes in their 
dust profiles. Indeed, if we do not exclude the bulges when applying linear fits to our dust 
profiles, the resulting dust scale-lengths are ~ 20% larger than the stellar ones for the whole 
sample, and ~ 30% larger if we focus on early-type spirals, whose bulges are more prominent. 

The bottom row in Fig. [8] shows the profiles of the dust light-to-mass ratio, which is 
essentially proportional to (U). Therefore, these profiles reflect radial changes in temperature 
which, in general, seems to decrease with radius. Interestingly, irregulars can reach larger 
dust luminosity surface densities than Sa-Sd spirals, in spite of having considerably lower 
dust mass densities, thanks to their dust being hotter (note their higher values of both Urain 
and 7). 

Note that we cannot speak of 'temperature' in the same way as when fitting FIR data 
with a double black-body SED, since the dust models used here consider the full temper- 
ature distribution for all grain sizes. Although the smallest grains can undergo extreme 
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temperature changes in small periods of time, large ones dominating the FIR emission can 
be considered to have reasonably steady temperatures T ^ 17f/^/^K (DL07). As a sanity 
check, for NGC 7331 we obtain a roughly flat profile with f/min ~ 2.5, implying T ^ 20K, 
the same value found by Thilker et al. (2007) when applying a double black-body fit. By 
analyzing the IR SEDs of different regions in M 81, Perez-Gonzalez et al. (2006) found that 
the temperature of the cold component varies from 22K in the central regions to ~ 15K 
at lOkpc from the center, which is consistent with our Umin profile ranging from 5 to 0.7 in 
the same relative radial range. 



4-3.3. Dust mass surface density and attenuation 

An initial comparison of the dust mass surface density profiles in Fig. [8] and the at- 
tenuation profiles in Fig. [1] shows that both quantities follow similar global trends along 
the Hubble sequence, peaking at Sb-Sbc spirals. However, the dust mass profiles show a 
smaller dispersion than the extinction profiles within each morphological-type bin (see also 
the histograms in Figs. [2] and [9]). 

To further analyze these differences, in Fig. [10] we plot the attenuation in the FUV 
as a function of the corresponding projected dust mass surface density (we cannot use the 
de-projected one, for the measured values of Apuv take into account all the dust along the 
line of sight). Note that since Apuv has been computed with the age-dependent calibration 
of Cortese et al. (2008), the varying relative contribution of young and old stars to the 
dust heating in galaxies with different SFHs is being taken into account. It could still be a 
problem in ellipticals and lenticulars, though. 

This plot reveals that, for a given value of the dust mass density, there is always a 
minimum associated attenuation. This value increases monotonically from almost zero for 
^M^dust proj ~ ^0 kpc~^ to roughly 3 magnitudes for S^/dust proj ~ ^0 kpc~^. This lower 
envelope for the distribution of data-points is delineated by Sc galaxies and later, while for 
the same amount of dust, earlier spirals present attenuations 1 or 2 mags higher. Using 
the boundary-fitting code described in Cardiel (2009), we obtained the following analytic 
parametrization for the minimum attenuation at a given dust column density: 

. r 5.16 - 2.57X + 0.327x2 if 4 < x < 7 

^"^^'-''^^l 0.1 if3<x<4 

where x = logS^f^^^^ in MQkpc"^. This fit is shown in Fig. [10] with a dashed line. 

The scatter in Apuv found for each given Sa/^^^^ can be attributed to the relative 
distribution of stars and dust within each galaxy. It is known that for geometries more com- 
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plex than a foreground homogeneous screen of dust, the distribution of dust grains relative 
to the heating sources overcomes the net amount of dust as the main factor determining the 
attenuation (Calzetti 2001). 

To illustrate this issue, in Fig. [11] we plot Apuv as a function of the far-infrared color. 
Data-points are clearly arranged in a wedge-like distribution, completely analogous to the one 
found by Dale et al. (2007) using the integrated photometry of the SINGS galaxies. These 
authors came to the conclusion that the location of galaxies in this diagram was closely 
related with their morphology in the 24 /im band. Galaxies with clumpy morphologies were 
found to have warm far-IR colors but low attenuations: while there are significant amounts 
of hot dust concentrated around heating sources in star-forming regions, many UV photons 
can find their way out of the galaxy through lines of sight clean of dust, thus decreasing 
the overall TIR-to-FUV ratio. In our radial analysis, Sc spirals and later clearly lie in the 
lower- right corner of the plot, thus being consistent with having a clumpy distribution of 
dust relative to stars. The presence of these kind of porous star-forming regions has been 
identified in nearby galaxies (Roussel et al. 2005), and their contribution to the TIR-to-FUV 
ratio may become progressively dominant with increasing galactocentric distance. 

In contrast. Dale et al. (2007) found that the upper- right region of the plot is populated 
by galaxies showing centrally concentrated emission at 24 /im. In our radially resolved anal- 
ysis the central regions of spirals, especially in Sb-Sbc ones, lie in this area of the diagram. 
This indicates that UV light has significant difficulties to avoid dust, thus leading to much 
larger extinction values. Finally, regions having cooler far-IR colors tend to exhibit interme- 
diate attenuations. This area of the plot was identified by Dale et al. (2007) as corresponding 
to galaxies with a smooth 24 /im morphology. This is also the case for our radial profiles, 
since the disks of Sbc spirals and earlier populate this zone of the plot. While not being as 
heavily extincted as in the central regions, UV photons in these disks still find it difficult 
to leak through dust, in comparison with clumpy Sc-Sd spirals, in which the attenuation is 
lower despite having similar dust mass surface densities. 

4.3.4- PAH abundance 

Over the years, a body of observational evidence has been collected showing that PAH 
emission is weak or even suppressed in low-metallicity systems (see e.g. Boselli et al. 1998; 
Sturm et al. 2000; Madden 2000; Madden et al. 2006; Hunt et al. 2005; Hogg et al. 2005; 
D07; Engelbracht et al. 2005, 2008). However, the physical mechanisms driving the observed 
correlation between PAH abundance and metallicity remain a subject of investigation. 
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Several authors have invoked an evolutionary origin for this trend (see e.g. Galliano et 
al. 2008; Dwek et al. 2009). PAHs are thought to form in the envelopes around carbon- 
rich AGB stars, being later injected into the ISM through a recycling process that takes 
place over timescales of a few Gyr. The remaining dust species are believed to condense 
out of ejecta from supernovae in much shorter timescales. Thus, most of the dust in young, 
low-metalhcity systems would be made out of SN-condensed species, whereas PAH particles 
would still require more time to form. Hence, the different speed of both processes would 
naturally lead to a progressive increase in (Jpah with time, explaining the observed PAH- 
metallicity trend. This scenario, however, should be considered cautiously, as many aspects 
of chemical evolution and dust formation still remain unclear. For instance, despite their low 
metallicities, some dwarf galaxies may harbor underlying stellar populations a few Gyrs old. 
Also, carbon dust might be also produced in the ejecta of supernovae, and its abundance 
relative to carbon dust due to AGB stars should be also taken into account. 

Conversely, other possible explanations rely on destructive mechanisms of PAH particles. 
Systems with younger stellar populations have harder UV radiation fields, given their larger 
fraction of ionizing OB stars. Furthermore, PAHs are likely to be less shielded by dust from 
these hard photons than in more metal-rich systems, since low-metallicity galaxies usually 
exhibit lower dust attenuation. There is observational evidence showing a progressive paucity 
of aromatic features as the hardness of the radiation field increases (Engelbracht et al. 2005, 
2008; Madden et al. 2006; Bendo et al. 2006, 2008; Smith et al. 2007; Gordon et al. 2008). 

It is not straightforward to disentangle the relative contribution of formation and de- 
struction effects to the observed PAH abundances. As we commented above, predominantly 
young systems which are expected to show a paucity of PAHs because of chemical evolution, 
also exhibit hard radiation fields that could destroy the already present ones. Besides, varia- 
tions in the aromatic emission could be due to changes in the ionization state of PAHs rather 
than to a real decrease in terms of abundance. Wu et al. (2006) found that the equivalent 
width of the aromatic features correlates better with a parameter simultaneously taking into 
account both the starlight hardness and the metallicity than with any of those parameters 
alone. 

To check whether the PAH-metallicity trend holds true in our radial analysis, in Fig. [T2] 
we plot the PAH abundance at each galactocentric distance as a function of the correspond- 
ing metallicity at the same radial bin, for the 22 galaxies in common with the sample of 
Moustakas et al. (2009, in preparation). As in Section [331 we only trust those metallicities 
interpolated within the spatial range covered by HII regions in the aforementioned study. As 
explained in Section 13. 3[ we adopt the metallicity values based on the Kobulnicky & Kewley 
(2004) calibration. Data points falling outside that range require uncertain extrapolations 
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of the metallicity gradients, and have been excluded from this analysis, as well as regions 
with excessive stellar emission at 8 fim (see Section U]2]). 

There is a general trend of qpAn to increase with the oxygen abundance for 12 + 
log(0/H) < 9. Note that in Sdm and irregular galaxies, the PAH abundance ranges from 
2.5% to less than 1%. However, there are not enough individual measurements of HIT re- 
gions in these galaxies to derive metallicity gradients, which precludes making use of those 
data-points for this plot. The available abundance data indicate that these galaxies have 
12 -|- log(0/H) < 8.6 (Moustakas et al. 2009, in preparation), so the monotonic trend be- 
tween gpAH and (0/H) continues down to lower abundances. This general trend is consistent 
with that found by D07 in the SINGS galaxies using galaxy-averaged PAH abundances and 
metallicities. The dashed line, which corresponds to a bisector linear fit applied to the data 
with 12 + log (0/H) < 9, shows that gpAH «: (0/H)°'®''^°'^ in this range of metallicities. 
The slope may vary within individual galaxies, and it could be different for lower oxygen 
abundances. We call attention to the fact that, up to now, integrated photometry usually 
required comparing normal spirals with dwarf galaxies in order to detect this trend, while 
radial profiles make it possible to see it within a more limited range of galaxy types because 
of the reduced scatter. 

At larger metallicities, however, the slope flattens and even reverses. In principle, the 
upper limit of the DL07 models for the PAH abundance (4.6%) could be the culprit, but we 
call attention to the fact that in most cases this flattening arises before that limit is reached, 
and it remains visible when recomputing the PAH abundances with the empirical recipes 
presented in Appendix IA.2I Interestingly, some models describing the stellar evolutionary 
effects on the chemical composition of dust predict that such a change in slope might actually 
occur (Galliano et al. 2008), and attribute it to the different carbon content of AGB stars 
with different masses —and hence with different lifetimes. Stars with masses between 1.2 and 
3.0 M0 are thought to end up as carbon-rich post-AGB stars, unlike low-mass stars, which 
remain oxygen-rich (see e.g. Garcia-Lario 2006 and references therein). Since the latter have 
longer lifetimes, they contribute later to the enrichment of the ISM. Moreover, the C/0 ratio 
in AGB stars also seems to depend on the stellar metallicity, with more metal-rich AGBs 
yielding lower amounts of carbon (see Dwek et al. 2009 and references therein). As a result, 
these models predict that the formation of carbon dust from AGB stars slows down with 
time, while the absolute abundance of SN-condensed dust species keeps rising. 

There are four galaxies in Fig. [12] that deserve special attention, for they have low 
PAH abundances compared to other galaxies with similar metallicities; moreover, the trend 
between those quantities is actually reversed within each one of these objects. Interestingly, 
they are the earliest spirals in the diagram (Sb and earlier). Therefore, they have presumably 
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undergone rapidly declining SFHs, with typical timescales of just a few Gyr (Gavazzi et al. 
2002), in contrast to late- type spirals, which exhibit longer timescales. If SFH is driving 
the progressive change of slope due to the different carbon content of AGBs with different 
masses and/or metallicities, this effect could be enhanced in systems with rapidly declining 
SFHs. 

In Fig. [13] we show how the PAH abundance changes with the observed and extinction- 
corrected (FUV — 3.6 /xm) color. The color of the radiation field illuminating the PAHs 
probably lies between the observed and extinction-corrected values because it may be par- 
tially reddened by dust. Regions with the reddest colors, usually the bulges of early-type 
spirals, exhibit low or intermediate PAH abundances, as we already saw in Fig. [81 But for 
(FUV — 3.6 yum)corr colors between and 4 mag (the range in which star-forming disks are 
typically found, as shown in Fig. [3]), PAHs progressively become less abundant in the outer 
and bluer regions, although the trend is noisy and mainly defined by a lower envelope. If 
photodestruction of PAH particles by intense UV fields overcame PAH formation by AGB 
stars, it would be conceivable to observe systematic differences in gpAH by sorting the data- 
points according to their FUV surface brightness (/xpuv) for a fixed (FUV — 3.6 fim) color, 
despite its limitations (see below). UV bright regions should populate the lower-left area of 
Fig. [13] if PAHs are being selectively destroyed in them, but they seem to have comparable 
PAH abundances as the remaining regions, since these brightest zones are usually the central, 
more evolved ones. 

Since we are lacking spatially resolved information from ionization tracers (such as the 
[NeIII]/[NeII] ratio), we are relying on the FUV surface brightness as a rough proxy for the 
level of ionization. This choice would be more justified in a pixel-to-pixel analysis, since /ipuv 
would then allow us to identify regions of massive star formation, but this analysis is beyond 
the scope of this paper. However, the usefulness of /ipuv in radial profiles as a ionization 
tracer is much more limited. It tends to select the central and brightest zones of the disks, 
where the evolutionary origin reflected by metallicity seems to be more important. In the 
outer regions, however, the azimuthally averaged FUV surface brightness is fainter, even if 
localized massive star formation is present. In other words, the declining FUV profiles do not 
imply that the UV field within individual star-forming regions becomes softer as we move 
away from the center of the galaxy. Indeed, it is worth noting that in M 33 the hardness of 
the UV radiation field inside single HII regions is seen to increase with galactocentric distance 
(Rubin et al. 2008). Since our PAH abundances are not limited to individual HII regions, 
but take into account all the aromatic emission within elliptical annuli, the corresponding 
ionization tracer should reflect the average radiation field heating all that dust, which is 
likely to be softer than the one present in the core of star-forming regions. With the available 
data, /ipuv is possibly our best choice, but in any case the results concerning the relative 
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importance of PAH destruction are not conclusive and should be treated with caution. 

The whole analysis presented in this section supports the idea that, at least in normal 
spiral galaxies and in spatial scales larger than those of single HIT regions, the observed 
trends between the PAH mass fraction and metallicity or (FUV — 3.6 /im) color may have an 
evolutionary origin. We emphasize, however, that while the varied timescales in which differ- 
ent elements are injected in the ISM may be partly responsible of the observed dependence 
of the PAH mass fraction on 0/H, other processes governing the balance of dust formation 
and destruction may be correlated with metallicity as well. 

In their analysis of individual HII regions in M 101, Gordon et al. (2008) find that the 
equivalent width of the aromatic features is better correlated with the level of ionization than 
with metallicity, and this seems to hold for the starburst galaxies studied by Engelbracht et 
al. (2008). Rather than representing a contradiction with other studies (including ours), these 
results reflect diverse behaviors associated with the different regimes of physical conditions 
being probed. While photodestruction of PAH particles by hard UV photons may overcome 
evolutionary effects in the core of HII regions and in starburst galaxies, it does not seem to 
play such an important role at larger spatial scales in our galaxies. Indeed, Gordon et al. 
(2008) note that this trend with the hardness of the radiation field only appears above a 
certain level of ionization. Below this level the trend seems to vanish, implying that PAHs 
are not being destroyed faster than the remaining dust species. This is likely the case in 
our azimuthally averaged profiles, where we are measuring aromatic features emerging not 
only from the very inside of HII regions, but also from the diffuse ISM, where PAHs are 
not necessarily exposed to such levels of ionization. From the SED decomposition shown in 
Fig. [6] we can see that, in general, the dust emission at 5.8 and 8.0 //m is dominated by the 
diffuse component. Therefore, at these spatial scales chemical evolution seems to act as the 
main driver of the observed trends, rather than selective PAH destruction. Recent analysis 
of the dust emission in NGC 0300 (Helou et al. 2004), NGC 4631 (Bendo et al. 2006) and a 
subset of 15 SINGS spirals (Bendo et al. 2008) seem to further support this general picture. 
As in Gordon et al. (2008), these authors found that the PAH emission at 8 //m and the hot 
dust radiation at 24 /im are correlated in scales of a few kpc, but this correlation breaks down 
within individual star-forming regions, where the PAH emission is less centrally peaked than 
the one contributed by very hot dust. 



4-3.5. Photometric 8 fim equivalent width 



Another important issue to address is how the PAH abundance —as a fraction of the 
total dust mass— relates to the equivalent width of the observed spectral aromatic features. 
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Although 2D IRS spectral maps are available for most of the SINGS galaxies, they do not 
usually cover the full spatial extent of these objects. In this regard, Engelbracht et al. (2008) 
propose a photometric method to estimate the equivalent width of the 7.7 /im feature using 
the stellar-subtracted flux at 8.0 /xm band, and a logarithmic interpolation of the stellar- 
subtracted fluxes at 4.5 and 24 /xm to estimate the underlying continuum. Note that, as 
pointed out by these authors, the 8 /im equivalent width (EW, either the spectroscopic or 
its photometric proxy) is a measurement of the abundance of PAHs relative to small grains 
or very hot larger ones, which contribute to the underlying continuum at those wavelengths. 
Translating the 8 //m EW into an abundance relative to the total dust mass can be challenging 
without a prior knowledge of the temperature distribution of the dust grains, as we will see 
below. 

It is not easy to obtain this photometric estimator for the 8 fim EW directly from our 
data. Although recovering the dust-only emission at 8.0 and 24 /im by means of the stellar 
factors presented in Section IT2] does not pose any problem, computing the non-stellar flux at 
4.5 /zm is not straightforward, since the emission at that band is almost entirely dominated by 
stars (the sample studied by Engelbracht et al. (2008) consists of starburst galaxies usually 
exhibiting less stellar contamination at that band). Indeed, in our sample the stellar emission 
at 4.5 nm, extrapolated from the flux at 3.6 /xm, is sometimes slightly larger than the actual 
observed values (partly due to photometric uncertainties, and to the fact that these stellar 
factors are average values). Alternatively, here we use the non-stellar fluxes corresponding 
to the best-fitting model for each galactocentric distance, instead of recovering them from 
the observed data. 

In Fig. [H] we plot the PAH abundance derived from the models against the photometric 
equivalent width estimated following the prescriptions given by Engelbracht et al. (2008). 
Both quantities are positively correlated, but there is a clear dependence on the fraction of 
dust exposed to intense starlight: at fixed gpAH, the photometric proxy for the equivalent 
width decreases as the fraction 7 of hot dust increases. This is most likely due to the effect 
that 7 has on the flux at 24 fim: as seen in previous sections, exposing a small amount of dust 
to very intense radiation fields can significantly increase the observed emission at that band, 
and thus increase the interpolated continuum at 8 /im. Therefore, the resulting equivalent 
width will result in smaller values. In conclusion, the 8 /xm EW may be used to compare the 
PAH mass fraction in samples of galaxies not presenting large variations in their FIR SEDs, 
but should be used with caution if this is not the case. 
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4- .3. 6. Dust-to-gas ratio 

The dust-to-gas ratio is an important quantity when studying the chemical enrichment 
of the ISM, since it tells us about the amount of metals that get locked up in the dust 
through the stellar yields. Indeed, a correlation has been widely shown between the dust- 
to-gas ratio and the gas-phase oxygen abundance (e.g. Lisenfeld & Ferrara 1998; Edmunds 
2001; James et al. 2002; Hirashita et al. 2002; D07). In their analysis of a small sample of 8 
nearby galaxies (including the MW), Issa et al. (1990) concluded that the dust-to-gas ratio 
(as traced by the Ay/Nn ratio) decreases with radius, following a trend with metallicity. 
A similar result was found by Boissier et al. (2004), who confirmed that a similar trend 
arises when the dust mass is directly estimated from the FIR emission, rather than from the 
attenuation. 

In order to obtain total gas profiles, we combined HI data from the THINGS survey 
with CO profiles compiled from the literature (see Table [1]) , which were transformed into 
H2 profiles by means of the metallicity-dependent CO-to-H2 conversion factor of Boselli et 
al. (2002). While this approach certainly constitutes an improvement over using a constant 
conversion factor, the derived H2 profiles should be still treated cautiously. The CO-to-H2 
conversion factor may depend on parameters like the cloud density and the excitation tem- 
perature of the CO, both of which could poentially vary systematically with galactocentric 
radius. Also, note that whereas the HI profiles were obtained in a consistent way with respect 
to the UV and IR data by measuring on the HI intensity maps degraded to the 160 /im res- 
olution, the CO profiles were directly interpolated at different radii from published profiles. 
Hence, the CO profiles might not necessarily match exactly the same spatial regions. Note 
also that since no metallicity gradients are available for NGC 4826 (Sab) and NGC 3627 (Sb), 
we used the CO-to-H2 conversion factor corresponding to the global if-band luminosity of 
each object (Boselli et al. 2002). 

In Fig.[T5]we show the radial dust-to-gas profiles for the SINGS-THINGS galaxies, sorted 
by morphological type. Those galaxies for which only HI data are available are plotted with 
a dashed-line, and their dust-to-gas ratios are therefore just upper limits, at least in the 
innermost regions where molecular gas is more likely to be found. 

In most objects the dust-to-gas ratio decreases with radius. In the central regions of Sb- 
Sd spirals it ranges between 10~^ and lO~^M0kpc~^, and it decreases faster in Sc-Sd ones, 
reaching values close to 10~^ kpc~^ or lower at the optical radius. NGC 3627 constitutes a 
notable exception, for its dust-to-gas ratio increases with radius. It is an interacting Sb spiral 
in the Leo triplet that gets redder in the outer regions, possibly because of SF quenching as a 
result of gas removal or exhaustion. NGC 2976 is a Sc galaxy whose dust to gas ratio remains 
roughly constant across the entire optical size. This quite unusual object exhibits two bright 
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knots of intense star formation located at almost symmetrical positions with respect to the 
center of the galaxy. 

In Fig. [TBk we plot the dust-to-gas ratio against the gas-phase oxygen abundance, for all 
the SINGS galaxies with available spatially resolved metallicities, HI and CO data. Points 
belonging to the same object have been connected for clarity. If we were to assume that the 
abundances of all heavy elements were proportional to the oxygen abundance and that all 
heavy elements condensed to form dust in the same way as in the MW, then the dust-to-gas 
ratio should scale proportionally to the oxygen abundance (D07): 

Mdust _ 0.01 (0/H) 

Mgas ^ 1.36 (0/H)mw ^ ' 

where 0.01 is the dust-to-hydrogen ratio of the MW and the factor 1.36 accounts for helium 
and heavier elements. The dashed line in Fig. [T6k marks this relation and shows that, at 
least to first order and for 12-|-log(0/H) > 8.9, the radial correlation between the dust-to-gas 
ratio and the metallicity follows this simple scaling law within a factor of ~ 2. In Fig. [T6b 
we plot the galaxies for which only HI data are available. Their values of Mdust/A^gas might 
be overestimated if significant amounts of molecular gas are present (which is likely the case 
in the central regions). 

We call attention to the fact that at low oxygen abundances the observed trend in panel 
(a) seems to be steeper than the one predicted by Eq. [6l A linear fit to the whole data-set 
yields: 

log(Mdust/Mg,3) = 5.63 + 2.45 x log(0/H) (7) 

In their analysis of the integrated photometry of the SINGS galaxies, D07 found that 
galaxies with and without available submillimeter data followed slightly different trends in 
the Mdnst/Mgs^ vs. metallicity plots, with the slope being a bit steeper when the dust masses 
had been obtained without submillimeter data. 

It is not clear, therefore, if the apparent departure from the MW scaling law seen in 
Fig.Hniat the lowest metallicities is real, or due to our cool dust masses being underestimated 
(the metallicity gradients might be also uncertain). Should this trend be real, it could link 
the behavior of the outermost, lower metal- abundant regions of spirals with the properties of 
dwarf galaxies, which exhibit lower dust-to-gas ratios than those predicted by extrapolating 
Eq. O Detailed models describing dust evolution in dwarf galaxies (see Edmunds 2001; 
Hirashita et al. 2002) usually invoke several mechanisms to either delay dust injection in 
the ISM or to eliminate it, via outflows caused by collective SN-driven winds. However, 
the extremely low dust-to-gas ratios of many dwarf galaxies could be partly due to their 
extended HI envelopes, where no dust emission is apparently seen (see e.g. Walter et al. 
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2007). Indeed, when Mdust/^gas is computed inside regions detected both in FIR and HI, 
the values derived are more consistent with Eq. [6] (DOT). 

The extent to which these explanations proposed for dwarf galaxies are also valid for the 
outer regions of spirals remains unclear. The radial decrease of the star-formation efficiency 
found in nearby spirals (Thilker et al. 2007; Leroy et al. 2008) could also partly explain the 
low dust-to-gas ratios seen in Fig. [161 The presence of large reservoirs of gas which has not 
been yet transformed into stars (and hence into dust) would drop the dust-to-gas ratio in 
the outer regions of these galaxies. 



5. Summary and conclusions 

We have measured multi-wavelength surface brightness profiles for the galaxies included 
in the SINGS sample, using UV, IR and HI data from GALEX, Spitzer (from SINGS) and 
VLA (from the THINGS survey), complemented with published CO data and metallicity 
gradients. The analysis of this panchromatic, spatially resolved data-set reveals important 
results regarding the dust properties of nearby galaxies: 

1. We have derived radial extinction profiles by means of a SFH-dependent calibration of 
the TIR-to-FUV ratio. In most galaxies the attenuation decreases with the galacto- 
centric distance. It is largest in Sb-Sbc galaxies, ranging from Afuv ~ 2.5 mag in the 
central regions to ~ 1.5 mag at the optical radius R25. As we move towards earlier 
or later Hubble types, the attenuation decreases and its radial variation becomes less 
pronounced. SO/a-Sab spirals typically have Apuv ~ 1.5 mag, and in Sdm and irregu- 
lar galaxies the extinction drops below 0.5 mag. However, the dispersion is rather large 
(~ 1 — 2 mag) within each bin of morphological types. This is the result of the different 
spatial distribution of stars and dust (see below). 

2. The TIR-to-FUV excess is correlated with the UV spectral slope (3, following an IRX-/3 
relation shifted with respect to the one found in starburst galaxies. The main driver 
of this departure seems to be the SFH, with more quiescent regions having redder 
UV slopes for a given TIR-to-FUV ratio. Indeed, the intrinsic (FUV — 3.6/im) color 
—a metric tracer of the current-to-past SFR— correlates well with the distance from 
the relation for starbursts. Our radial analysis reveals that this shift is particularly 
evident in bulges, as well as in the outer regions of anemic spirals with clear signs 
of star-formation quenching. However, second-order factors might also play a role, 
since the dependence on the SFH gets completely blurred unless a wide range of SFHs 
is explored. We also provide a calibration to estimate the UV attenuation from the 
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UV slope when FIR data are unavailable, which is statistically valid for 'normal' star- 
forming galaxies. 

3. By combining our data with oxygen abundance gradients, we have analyzed the influ- 
ence of metallicity on attenuation. Both quantities are clearly correlated, with metal- 
richer regions exhibiting larger attenuations. Inclination plays an important role here; 
at fixed metallicity, the extinction is larger in galaxies that are closer to edge-on. The 
same happens with the UV color: more inclined galaxies appear redder than nearly 
face-on ones. Both effects are roughly balanced when combining both quantities in the 
IRX-/3 plot, so that data-points are just shifted along the IRX-/3 relationship. 

4. We have fit the FIR profiles with the dust models of Draine & Li (2007), deriving the 
radial distribution of the PAH abundance, the dust mass and luminosity surface den- 
sities and the properties of the heating starlight. The fraction of the dust mass in the 
form of PAHs (^pah) varies widely along the Hubble sequence, reaching a maximum of 
~ 3 — 4.5% in Sb-Sbc spirals, and almost disappearing (< 1%) in the outer regions of 
Sdm and irregular galaxies. Indeed, the PAH abundance typically decreases with radius 
in Sc spirals and later; in earlier types, however, the profiles tend to flatten and even 
reverse. The PAH abundance increases with metaUicity for 8.5 < 12 -|- log(0/H) < 9, 
and the trend progressively flattens and eventually reverses at larger oxygen abun- 
dance. This is in qualitative agreement with the predictions of models coupling stellar 
evolution with dust formation. Although this dependence on metallicity has already 
been found when comparing normal spirals with dwarf galaxies, using radial profiles 
leads to a considerable reduction of the dispersion and carrying out the study in a 
narrower range in metallicity, without the need to include low-met allicity dwarfs. 

5. The average intensity of the heating starlight decreases monotonically with the galac- 
tocentric distance, and in the outer regions it usually reaches similar values as the local 
MW radiation field, corresponding to cold dust temperatures of ~ 20 K. In order to 
explain the observed 24 fim fiux densities, it is required that about 1% of the total dust 
mass must be exposed to very intense radiation fields. Such a small amount of dust 
can account for ~ 10% of the total dust IR luminosity. 

6. The dust mass surface density is largest in Sb-Sbc spirals, varying from 10^ to 10^ Mq kpc 
from the center to the R25 radius. Dust is less abundant in earlier spirals, typically 
having surface densities below 10^ Mq kpc~^, and also in Sdm and later, where it drops 
below lO^Mgkpc"^. The dust profiles are exponential, and have a radial scale-length 
that holds remarkably constant (0.2-0.3 R25) from Sb to Sd galaxies. The dust pro- 
files in so/a and Sab spirals usually present a central depletion, but have similar radial 
scale-lengths than Sb-Sd ones in the outer exponential regions. Compared to the radial 
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distribution of old stars, as traced by the 3.6 //m profiles, the dust radial scale-length 
is, on average, ~ 10% larger than the stellar one. 

7. The relative spatial distribution of stars and dust plays a crucial role in determining 
the attenuation. There is a wide range of attenuations (1-2 mag) that can be found 
for a given observed (i.e. projected) dust mass surface density. The minimum values 
of attenuation are always found in Sc galaxies and later, while for the same amount 
of dust earlier spirals (especially Sb-Sbc ones) exhibit larger extinctions. In late-type 
spirals, dust and stars are arranged in a clumpy geometry: while some heating sources 
can be heavily attenuated, many UV photons can leak through lines of sight clean 
of dust, thus decreasing the overall attenuation. This does not happen so easily in 
the disks of Sbc spirals and earlier, where the same dust surface density yields larger 
attenuations. 

8. By merging the dust mass profiles with HI data from the THINGS survey and CO data 
from the literature, we have studied the radial variation of the dust-to-gas ratio, which 
is seen to decrease by an order of magnitude from the center to the edge of the optical 
disk of each galaxy. Typical values range from 0.1 - 0.01 in the central regions to 0.01 - 
0.001 in the outer ones, with early-type spirals having larger amounts of dust than 
late- type ones with a similar gas content. The dust-to-gas ratio is clearly correlated 
with the gas-phase oxygen abundance, being larger in the central, metal-rich zones of 
the disks than in the outer ones. To first order, for 12 + log(0/H > 8.9) the relation 
can be described with a simple scaling law: M^ust/^gas (0/H), which would imply 
that most condensable elements tend to form dust in a roughly similar way in spiral 
galaxies of different morphological types. However, at lower metallicities the dust-to- 
gas ratios are systematically below this simple relation. This could imply that large 
amounts of gas that has not yet undergone star-formation activity reside in the outer 
regions of spiral disks, although the presence of dust cooler than 15 K could increase 
our dust-to-gas ratios. 
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A. Empirical estimations of the model parameters 

Fitting the infrared SEDs of star-forming galaxies with detailed models for the dust 
emission constitutes the proper way to analyze the physical properties of the interstellar 
dust. However, it is also desirable to have empirical calibrations from which one can roughly 
estimate these properties directly from observed quantities. Here we provide some basic 
equations for that purpose. In Table [5] we provide the residual rms of each one of these 
recipes. 



A.l. Estimating 7 

Draine & Li (2007) already provide a useful fit to estimate /(L^; U > 10^), that is, the 
fraction of the total dust luminosity contributed by high- intensity regions. They define the 
following quantity: 

P.-0,14a„. <7^"-°-"<rf>- (Al) 

where the nonstellar flux densities F^^ can be estimated using the stellar emission corrections 
given in Section 14. 2[ The numerator is the 24 /xm flux coming from large grains heated by 
very intense radiation fields. The factor 0.14(i/F"'')8,o accounts for the emission at 24 /im 
contributed by stochastically heated PAHs, since this emission is not necessarily linked to 
high-intensity regions. DL07 showed that /(L^; U > 10^) can be estimated using the follow- 
ing fitting function: 

fiU, U > 102)(%) = 105(P24 - 0.14P8.0 - 0.035)°-^^ (A2) 



While DL07 showed that this function is certainly a good proxy for the fraction of the 
dust luminosity coming from high-intensity regions. Fig. [TTk shows that this is also the case 
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for our radial profiles. The dispersion in the lower region is due to differences in ?7min and 

QPAH- 

The corresponding mass fraction enclosed in these regions is also tightly correlated with 
P24 — 0.14P8 0, as can be seen in panel (b). Following DL07, we performed a similar fit to 
obtain an estimator for 7: 

7 (%) = 46.2(P24 - 0.14P8.0 - 0.023)^-2^ (A3) 

A. 2. Estimating gpAH 

The emission at 8.0 fim depends on both the PAH abundance and the amount of dust 
exposed to high- intensity starlight. Therefore, one would expect that gpAH could be obtained 
from Pg.o and P24 — 0.14P8.0, the latter accounting for differences in 7, as seen above. This is 
confirmed in Fig. [T8k . where we can see that the PAH abundance correlates well with Pg.o, 
with a second order dependence on P24 — 0.14P8.0. As expected, at fixed gpAH the emission 
at 8.0 fim —relative to the emission at longer wavelengths— increases with 7. 

A close inspection of data-points distribution reveals that a certain linear combination 
of Pg.o and P24 would allow estimating gpAH with a single fit. We computed this linear 
combination by means of a Principal Component Analysis. The results are shown in panel 
(b), where we can see that the PAH abundance can be readily estimated using the following 
fitting function: 

gpAH (%) = 13.03(P8.o - O.38P24 + 0.022)°-^^ (A4) 

Since the largest value of the PAH abundance derived from the models is gpAH = 4.6%, 
there are many data-points saturating that region of the plot —their actual PAH abundances 
being possibly larger— so we excluded those points in this analysis. For the same reason, 
points with gpAH = 0.4% were also excluded, since this is the lower limit of the dust models. 
Nevertheless, our fit should still allow us to recover PAH abundances a few percent larger 
than the current upper limit, since there is no reason a priori for this trend to break just 
above that value. 

A. 3. Estimating (U) 

The scale factor of the dust- weighted mean starlight intensity, {U), depends strongly 
on the shape of the FIR SED. In Fig. [19] we show that this parameter is linearly correlated 
—in logarithmic scale— with P70 = {vFy)^^/ {i'Fy)iQQ, which is sensitive to the temperature 
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of the largest grains dominating the FIR emission. A hnear fit yields the following empirical 
function: 



\og{U) = 0.468 + 1.801 log i?7o (A5) 

The starlight intensity of the diffuse component, f/min, is not so tightly correlated with 
i?7o, since the effect of 7 must be taken into account. Once {U) and 7 have been estimated 
using the formulae above, Umin can be computed from Eq. 33 in DL07. 

Note that, as explained in Section IH the smallest value of Umm considered in the model- 
fitting procedure is 0.7, as DL07 suggest to proceed when lacking submillimeter data. As 
a consequence, the resulting fits also have {U) > 0.7. For this reason, some data-points 
are clustered in the lower region of Fig. and have not been taken into account when 
performing the fit. 



A. 4. Estimating the dust mass 

We can estimate the dust mass from the mean starlight intensity (U) and the dust 
luminosity, since 

Ldnst = Po{U)Mdnst (A6) 

where Pq is the power per unit of mass radiated by the dust, when exposed to a radiation 
field equal to the local MW one (i.e. with U = 1). According to the DL07 models, differences 
in the PAH abundance can introduce variations of ~ 2% in this parameter. Nevertheless, a 
constant average value of Pq ~ I37L0/M0 can be safely adopted. 

The dust luminosity can be estimated directly from the photometric data using the 
weighted sum proposed by DL07. By substituting that sum and the empirical fit for (U) 
(Eq. lASp into Eq. IA6[ we can derive the dust mass from the IRAC and MIPS fiux densities 
and the distance to the source: 

Mdust ^ 47r f D y f {uF,)7o y'-''' 

Mq 1.616 X 10-13 [upcj \{iyF,)ieoJ 

0.95(z/F;^)8.o + 1.150(z/F-)24 + (^/P.)70 + (^^i^.)i60 , . 

^ ] 5 l^'J 

erg s~^ cm~"' 

As in previous sections, the non-stellar fiux densities F"'^ can be obtained from the 
observed ones by extrapolating the stellar continuum from the observed 3.6 /xm fiux density 
(see Section \A.2\i . To check the validity of this estimator, we have applied it to our radial 
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profiles. In Fig. [20b we compare the dust mass surface densities obtained from the model- 
fitting procedure with the ones estimated through Eq. IA7[ The agreement is excellent, with 
a scatter of only 9%. The most deviant points are those with the smallest value of the diffuse 
starlight intensity considered throughout this work (see Appendix IA.3p . Note that there is 
a small offset between the dust masses derived from Eq. IA7I and the actual ones obtained 
from the model fitting, due to the fact that the Ldust estimator slightly overestimates the 
actual Ldust from the model for the particular range of values of {U) and gpAH in our galaxies. 
Although Eq. IA7I is intended to be more general, if one wishes to account for this fact, the 
dust masses derived from Eq. IA7I should be multiplied by 0.95 (see Table [5]). 

In the absence of IRAC data, we can still compute the dust mass simply by substituting 
the Ltir estimator of DL07 (i.e. the numerator of the last term in Eq. IA7p with the MIPS- 
only calibration of Dale & Helou (2002). We can also skip the stellar emission correction for 
the 24 /im band, since it is almost negligible. Hence, we have: 



In Fig. [20b we see that this expression also constitutes a good metric tracer for the dust 
mass. Note, however, that the dust masses derived in this paper have been obtained in the 
absence of spatially resolved submillimeter data, and therefore the amount of dust colder 
than 15 K might not be strongly constrained. When analyzing the integrated properties of 
the SINGS galaxies, D07 showed that omitting submillimeter data can lead to an uncertainty 
of ~ 50% in the dust mass, and this should be taken into account when using the previous 
equations to directly compute the dust masses. 



When using the DL07 dust models to fit our radial infrared SEDs, we have relied only 
on IRAC and MIPS data. Hence, it is worth analyzing the possible systematics effects that 
may arise from using a limited number of bands to constrain the models. D07 made use of 
a larger set of multiwavelength data to study the global properties of the SINGS galaxies. 
Besides IRAC and MIPS observations, they included in their fit data from 2MASS, IRAS, 
SCUBA (for 17 galaxies) and the IRS "blue peak-up" detector array onboard Spitzer (for 9 
galaxies). While useful for global analysis, the use of these additional data sets in a spatially 
resolved analysis poses some problems regarding resolution, depth and heterogeneous data 
quality. 




X 



1.559(;/F,)24 + 0.7686(;/F,)7o + 1.347(z/F,)i6o 
erg s~^ cm~2 



(A8) 



B. Possible systematic effects in the dust models 
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To assess the validity of our methodology, we have compared our values of the model 
parameters with the ones derived by DOT for each galaxy as a whole. In principle, given 
a radial profile for each parameter in a certain galaxy, it is not straightforward to reduce 
'a posteriori' the spatially resolved information into a single data-point representative of 
the whole galaxy. A more reliable and robust procedure consists of applying the empirical 
formulae presented in Appendix |X] to the IRAC and MIPS integrated photometry used in 
D07. These recipes are affected by the same systematic effects —if any— as the values 
directly yielded by our model fitting procedure, since they were explicitly calibrated using 
those data. Should there be any strong biases in our values, they should arise when comparing 
our estimates with the values computed by D07. 

Table quotes the offsets and scatter between our values and those published in D07. 
Note that part of the scatter is due to the intrinsic rms of each recipe. In Fig. [2T]we compare 
the D07 values with the ones estimated using our recipes. Panel (a) shows that our PAH 
abundances might be overestimated at very low abundances, and underestimated at the 
largest ones. The overall relative scatter is ~ 30%. 

In panel (b) we compare the values of the fraction of dust mass heated by very intense 
starlight. The correlation is excellent, although our values are systematically 6% lower, with 
a rms of ~ 25%. 

The average intensity of the radiation field heating the dust that we get from our SEDs 
is just 3% lower than the one from D07, with a scatter of 22%, as can be seen in panel 
(c). The most deviant points contributing to the dispersion are those with available SCUBA 
data, for which we get lower dust temperatures. 

As a result, our dust masses appear to be just slightly overestimated in comparison with 
those obtained by D07 [see panel (d)]. For galaxies lacking SCUBA data, our masses are 
only 6% larger, with a rms of ~ 22%. In galaxies with available submillimeter data, the 
offset increases up to 25%, with a scatter of 40%. 

In brief, the correlations shown in Fig. [5T] are good enough to assert that the results 
presented in this paper do not suffer from any strong bias due to our restricted wavelength 
range that could compromise the validity of our conclusions. 
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Table 1. Sample 





1^^.2000 


DEC2OOO 


2a 


2b 


P.A. 


E(B— V) 


T 


dist 




THINGS 


CO 








( !^ Vr'tTl 1 Tl 1 


( !^ vr'iTi 1 n 1 

I CLl. \ulLIIII 1 


(deg) 




type 












(2) 




(4) 


(5) 


(6) 


(7) 


(8) 


(9) 




cm 




NGC 0024 


00 09 56.5 


-24 57 47.3 


5.8 


1.3 


46 


0.020 


5 


8.2 


yes 


no 




NGC 0337 


00 59 50.1 


-07 34 40.7 


2.9 


1.8 


310 


0.112 


7 


25 


yes 


no 




NGC 0628 


01 36 41.8 


15 47 00.5 


10.5 


9.5 


25 


0.070 


5 


11 


yes 


yes 


1 


NGC 0855 


02 14 03.6 


27 52 37.8 


2.6 


1.0 


60 


0.072 


-5 


9.7 


yes 


no 




NGC 0925 


02 27 16.9 


33 34 45.0 


10.5 


5.9 


282 


0.076 


7 


9.3 


yes 


yes 


2 


NGC 1097 


02 46 19.1 


-30 16 29.7 


9.3 


6.3 


310 


0.027 


3 


15 


yes 


no 


3 


NGC 1291 


03 17 18.6 


-41 06 29.1 


9.8 


8.1 


345 


0.013 





9.7 


yes 


no 




NGC 1316 


03 22 41.7 


-37 12 29.6 


12.0 


8.5 


50 


0.021 


-2 


19 


yes 


no 




NGC 1512 


04 03 54.3 


-43 20 55.9 


8.9 


5.6 


90 


0.011 


1 


10 


yes 


no 




NGC 1566 


04 20 00.4 


-54 56 16.1 


8.3 


6.6 


60 


0.009 


4 


17 


yes 


no 


4 


NGC 1705 


04 54 13.5 


-53 21 39.8 


1.9 


1.4 


50 


0.008 


11 


5.1 


yes 


no 




NGC 2403 


07 36 51.4 


65 36 09.2 


21.9 


12.3 


307 


0.040 


6 


3.2 


yes 


yes 


3 


Holmbcrg II 


08 19 05.0 


70 43 12.1 


7.9 


6.3 


15 


0.032 


10 


3.4 


yes 


yes 




DDO 053 


08 34 07.2 


66 10 54.0 


1.5 


1.3 


300 


0.037 


10 


3.6 


yes 


yes 




NGC 2841 


09 22 02.6 


50 58 35.5 


8.1 


3.5 


327 


0.016 


3 


14 


yes 


yes 


5 


NGC 2915 


09 26 11.5 


-76 37 34.8 


1.9 


1.0 


309 


0.275 


90 


3.8 


yes 


no 




Holmberg I 


09 40 32.3 


71 10 56.0 


3.6 


3.0 


360 


0.048 


10 


3.8 


yes 


yes 




NGC 2976 


09 47 15.5 


67 54 59.0 


5.9 


2.7 


323 


0.069 


5 


3.6 


yes 


yes 




NGC 3049 


09 54 49.7 


09 16 17.9 


2.2 


1.4 


25 


0.038 


2 


22 


NUV 


no 




NGC 3031 


09 55 33.2 


69 03 55.1 


26.9 


14.1 


337 


0.080 


2 


3.6 


yes 


yes 


2 


NGC 3190 


10 18 05.6 


21 49 55.0 


4.4 


1.5 


305 


0.025 


1 


17 


yes 


no 




NGC 3184 


10 18 17.0 


41 25 28.0 


7.4 


6.9 


135 


0.017 


6 


8.6 


yes 


no 




NGC 3198 


10 19 54.9 


45 32 59.0 


8.5 


3.3 


35 


0.012 


5 


17 


yes 


yes 




IC 2574 


10 28 23.5 


68 24 43.7 


13.2 


5.4 


50 


0.036 


9 


4.0 


yes 


yes 




NGC 3351 


10 43 57.7 


11 42 13.0 


7.4 


5.0 


13 


0.028 


3 


12 


yes 


yes 




NGC 3521 


11 05 48.6 


-00 02 09.1 


11.0 


5.1 


343 


0.058 


4 


9.0 


yes 


yes 


1 


NGC 3621 


11 18 16.5 


-32 48 50.6 


12.3 


7.1 


339 


0.080 


7 


8.3 


yes 


yes 




NGC 3627 


11 20 15.0 


12 59 29.6 


9.1 


4.2 


353 


0.032 


3 


9.1 


yes 


yes 


1 


NGC 3938 


11 52 49.4 


44 07 14.6 


5.4 


4.9 


15 


0.021 


5 


12 


NUV 


no 




NGC 4125 


12 08 06.0 


65 10 26.9 


5.8 


3.2 


275 


0.019 


-5 


21 


NUV 


no 




NGC 4236 


12 16 42.1 


69 27 45.3 


21.9 


7.2 


342 


0.015 


8 


4.5 


yes 


no 




NGC 4254 


12 18 49.6 


14 24 59.4 


5.4 


4.7 


35 


0.039 


5 


17 


NUV 


no 




NGC 4321 


12 22 54.9 


15 49 20.6 


7.4 


6.3 


30 


0.026 


4 


18 


NUV 


no 





Table 1 — Continued 



Object name 


RA2000 


DEC200() 


2a 


2b 


P.A. 


E(B-V) 


T 


dist 


GALEX 


THINGS 


CO 




(h;m:s) 


(d:m:s) 


(arcmin) 


(arcmin) 


(deg) 


(mag) 


type 


(Mpc) 


data 


data 


data 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 




1 9 9Q on ft 


J- ( uo uo.o 


^ 9 
o.z 


o.y 


000 


n9R 

U.UZo 


2 


1 7 


NT TV 
IN U V 


no 






1 2 S4 27 1 


02 1 1 1 fS 4 

KJ ^ -Li J_0»^ 


7.6 


3.2 


310 


0.018 


4 


15 


yes 










97 ^^7 '^'^ 1 


10.7 


4.4 


330 


0.018 


6 


17 


yes 








19 "^R 4Q K 


1 HQ 4R 


9.5 


4.4 


23 


0.046 


2 


17 


yes 




1 6 


vjr w i y 


1 9 '^7 4"^ 


1 1 4Q 1 


o.y 


4 7 


97^^ 
z / 


n 041 


Q 




1 7 


yes 


no 






1 9 '^Q "^Q 4 


— 11 '^7 9^0 


8.7 


3.5 


90 


0.051 


1 


9.1 


yes 








1 9 41 "^9 7 


41 1 R 9^; 4 


2 2 


1 Q 

-L .y 


oou 


01 s 


Q 

y 


y .0 


yes 


no 






19 4 9 riR n 

IZ "iZ Uo.U 


QO QO 90 /l 

oz oz zy.^ 


10.0 


9 7 

Z. f 


SR 


017 
U.Ul 1 


7 


Q 

y .u 


yes 






NOP 479^ 


19 "^0 9R fi 

J- Z OVJ ZiJ.U 


9^^ "^n 09 7 
zo ou UZ. ( 


10.7 


7.6 


35 


0.012 


2 


17 


yes 








1 9 f^n '^'^ 1 
J- z ou 00. ± 


41 07 1 fi 


11.2 


9.1 


285 


0.018 


2 


5.2 


yes 


yes 




IN V_I ^O^XJ 


iz ou ^0.0 


91 40 i^l Q 
zx ox.y 


10 


4 
0.^ 


zyo 


041 


2 


1 7 


yes 


yes 




TvTPP t^fT^^ 
iN'orl_/ OUOO 


1 Q 1 Q 97 C 


oa Qc QQ n 
OD OO 00. U 


10 7 


O.U 


1 7n 


n ni 1 

U.Ul 1 





1 


yes 


no 




NGC 5055 


13 15 49.3 


42 01 45.4 


12.6 


7.2 


285 


0.018 


4 


8.2 


yes 


yes 


1 


NGC 5194 t 


13 29 52.7 


47 11 42.6 


11.2 


9.0 


360 


0.035 


4 


8.4 


yes 


yes 


7 


TOL 89 


14 01 21.6 


-33 03 49.6 


2.8 


1.7 


352 


0.066 


8.1 


16 


yes 


no 




NGC 5408 


14 03 20.9 


-41 22 40.0 


1.6 


0.8 


12 


0.069 


9.7 


4.5 


no 


no 




NGC 5474 


14 05 01.6 


53 39 44.0 


4.8 


4.3 


360 


0.011 


6 


6.8 


yes 


no 




NGC 5713 


14 40 11.5 


-00 17 21.2 


2.8 


2.5 


10 


0.039 


4 


27 


yes 


no 




NGC 5866 


15 06 29.6 


55 45 47.9 


4.7 


1.9 


308 


0.013 


-1 


15 


yes 


no 




IC 4710 


18 28 .38.0 


-66 58 56.0 


3.6 


2.8 


5 


0.089 


9 


8.5 


yes 


no 




NGC 6822 


19 44 56.6 


-14 47 21.4 


15.5 


13.5 


7 


0.236 


10 


0.60 


yes 


no 




NGC 6946 


20 34 52.3 


60 09 14.2 


11.5 


9.8 


75 


0.342 


6 


5.5 


yes 


no 




NGC 7331 


22 37 04.1 


34 24 56.3 


10.5 


3.7 


351 


0.091 


3 


15 


yes 


yes 


5 


NGC 7793 


23 57 49.8 


-32 35 27.7 


9.3 


6.3 


278 


0.019 


7 


2.0 


yes 


yes 





Note. — Sample. (1): Galaxy name. (2): RA(J2000) of the galaxy center. (3): DEC (J2000) of the galaxy center. (4), (5): Apparent 
major and minor isophotal diameters at /Us=25 mag arcsec"^ from the RC3 catalog. (6): Position angle from RC3. JThe PA and axis 
ratio for NGC 5194 differ from those in the RC3, which are affected by the presence of NGC 5195. (7): Galactic color excess from Schlegel 

et al. (1998). (8): Morphological type T as given in the RC3 catalog. (9): Distance to the galaxy, rounded to the nearest Mpc when larger 
than 10 Mpc, taken from Gil dc Paz et al. (2007) and Kennicutt et al. (2003). (10): Available GALEX images. (11); Available HI maps 
from THINGS. (12): References for the CO data. 1: Regan et al. (2001). 2: Sage (1993). 3: Young et al. (1995). 4: Bajaja et al. (1995). 
5: Young & Scoville (1982). 6: Kenney & Young (1988). 7: Paglione et al. (2001). 



Table 2. Degraded GALEX, IRAC and MIPS profiles. 



r f A'FUV Mnuv M3.6/im A'4.5/im At5.8/im /i3.6/im At24/im M70/im /iieo/im 

(arcsec) (kpc) (mag/D) (mag/D) (mag/D) (mag/D) (mag/D) (mag/D) (mag/D) (mag/D) (mag/D) 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 



NGC 7331 




48 
96 
144 
192 
240 
288 
336 
384 
432 



0.0 
3.5 
6.9 
10.4 
13.8 
17.3 
20.7 
24.2 
27.7 
31.1 



24.99±0.02 
25.08±0.01 
25.31±0.01 
25.62±0.01 
26.07±0.01 
26.52±0.01 
27.30±0.01 
28.14±0.03 
28.88±0.05 
29.43±0.07 



23.74±0.01 
24.01±0.00 
24.53±0.01 
24.97±0.01 
25.47±0.01 
25.99±0.02 
26.74±0.04 
27.51±0.08 
28.15±0.15 
28.58±0.22 



18.41±0.01 
18.99±0.01 
20.03±0.01 
20.78±0.01 
21.39±0.01 
22.06±0.02 
22.77±0.04 
23.32±0.07 
23.69±0.09 
23.98±0.12 



18.88±0.01 
19.45±0.01 
20.47±0.01 
21.22±0.01 
21.84±0.02 
22.50±0.03 
23.20±0.06 
23.75±0.10 
24.12±0.14 
24.42±0.18 



18.48±0.01 
18.90±0.01 
19.74±0.01 
20.52±0.02 
21.20±0.04 
21.88±0.07 
22.66±0.14 
23.33±0.26 
23.79±0.39 
24.15±0.54 



17.81±0.01 
18.12±0.01 
18.86±0.01 
19.65±0.01 
20.36±0.02 
21.07±0.04 
21.93±0.09 
22.68±0.17 
23.18±0.28 
23.53±0.38 



17.82±0.01 
18.04±0.01 
18.83±0.01 
19.78±0.01 
20.59±0.02 
21.34±0.04 
22.28±0.10 
23.11±0.22 
23.69±0.37 
23.96±0.47 



15.08±0.01 
15.16±0.01 
15.71±0.01 
16.56±0.01 
17.39±0.03 
18.10±0.05 
18.79±0.09 
19.44±0.16 
19.81±0.22 
20.06±0.28 



14.11±0.01 
14.15±0.01 
14.82±0.01 
15.60±0.01 
16.28±0.01 
17.03±0.02 
17.82±0.04 
18.43±0.07 
18.95±0.22 
19.38±0.17 



Note. — Surface brightness profiles measured on the GALEX, IRAC and MIPS images, after matching their PSFs to that of 
the 160 micron band. Only the profiles for NGC 7331 are shown as an example in the printed version of the jornal. Check the 
on-line edition for the full table. (1): Radius along the semi-major axis in arcseconds. (2): Radius along the semimajor axis in kpc. 
(3)-(ll): Surface brightness at different wavelengths. All values are in AB mag arcsec"'^. The uncertainties include photometric and 
background errors, but not zero-point uncertainties, which are added in quadrature when needed. 
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Tablc 3. Extinction profiles 



B05 C08 

r r log(LTIR/iFUv) log(LTIR/i'NUv) ^FUV ^NUV ^FUV ^NUV 

(arcsec) (kpc) (mag) (mag) (mag) (mag) 

(1) (2) (3) (4) (5) (6) (7) (8) 



NGC 7331 






0.0 


1 


76 


± 


0.07 


1 


44 


± 


0.07 


3 52+" 11 


9 ..c+O.lO 


2 45+"-"' 
^•^3-0.10 


1 


-"^-0.10 


48 


3.5 


1 


76 


± 


0.07 


1 


50 


± 


0.07 


3.50«:il 




9 71 +0.11 
^- '1-0.11 


2 


91+0.11 
^-^-0.11 


96 


6.9 


1 


59 


± 


0.07 


1 


45 


± 


0.07 




9 qo+0.10 

^••^"-o.io 


2.691°;}} 


2 


9R + 0.12 
^"-0.12 


144 


10.4 


1 


39 


± 


0.07 


1 


30 


± 


0.07 


9 7C+0.10 
^- '°-0.10 


2 OQ+O O^ 


2.4liri[! 


2 


Q1+O.IO 
'^1-0.10 


192 


13.8 


1 


27 




0.07 


1 


21 


± 


0.07 


^•■^l-O.O!) 


1 01 +0-09 
l-yi-0.09 


9 91 +0.10 
^•^1-0.10 


1 


QO+0.10 

^■^-0.10 


240 


17.3 


1 


16 


± 


0.07 


1 


12 


± 


0.07 




1 7C+0.09 
1- "J-0.09 


2.04i»i;:5 


1 


pq+0.09 
"■'-0.09 


288 


20.7 


1 


16 




0.07 


1 


11 


± 


0.07 




1 74+0.09 
1- '^-0.09 


2-03+s!-};; 


1 


66+0-10 
"'^-0.10 


336 


24.2 


1 


23 


± 


0.07 


1 


15 


± 


0.08 


9 44+0.10 
^•4^-0.10 


-I 09+0.11 

l-°^-o.ii 


2 ii+''-i" 

^■11-0. 10 


1 


7Q+0.11 
"'-0.11 


384 


27.7 


1 


34 


± 


0.08 


1 


22 


± 


0.10 


^■"^-0.12 


1 94+"-i4 


9 90+0. 12 
^■^■J-0.12 


1 


70+0.14 

'^-0.14 


432 


31.1 


1 


43 


± 


0.09 


1 


26 


± 


0.12 


9 09+0. 14 


2 01+°-l* 


9 09+0.14 
''•"J^-0.14 


1 


00+0. 18 

"•^-c.is 



Note. — Radial extinction profiles. (1): Radius along semimajor axis in arcsec. (2): Radius along semimajor 
axis in kpc. (3): TIR-to-FUV ratio. (4): TIR-to-NUV ratio. (5), (6): Extinction at FUV and NUV computed 
with the SFH-indcpcndcnt fits of Buat et al. (2005). (7), (8): Extinction at FUV and NUV computed with the 
SFH-dependent prescriptions of Cortese et al. (2008). 
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Table 4. Model parameters profiles 



r- T (?PAH 7 

(arcsec) (kpc) (%) (%) (%) (Lq kpc'^) (Mokpc-^) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 



NGC 7331 
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-0.3 
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Note. — Radial profiles of the model parameters. (1): Radius along semimajor axis in axcsec. (2): Radius along scmimajor 
axis in kpc. (3): Fraction of the dust mass in the form of PAHs. This value is not reliable in regions marked with a dagger (f), 
where the dust contribution to the observed 8/im flux is lees than half of the stellar emission at that band. (4): Fraction of 
the dust mass heated by very intense starlight. (5): Fra<;tion of the dust luminosity contributed by regions heated by starlight 
■with U > 10-^. (6): Minimum value for the starlight scale factor, in units of the local MW radiation field. (7): Dust- weighted 
average scale factor for the starlight intensity. (8): Dust luminosity surface density, corrected for inclination. (9): Dust mass 
surface density, corrected for inclination. (10): Dust-to-gas ratio. The dust mass is computed as Mgas = 1.36 X (Mhi -|- M^^), 
or Mgas = 1.36 X Mm when CO data are not available. In the latter case, the values are marked with a double dagger (|), and 
should be treated as upper limits (at least in the inner regions). 
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Table 5. Accuracy of the recipes 



Parameter^ 


Empirical fit^ 




Comparison with D07 




no SCUBA^ 


SCUBA* 


offset 


rms 


ofTset 


rms 


offset 


rms 


<?PAH 


0% 


±11% 


+2.7% 


±30% 


-1.1% 


±14% 


7 


0% 


±25% 


-5.7% 


±23% 


-6.4% 


±27% 


(U) 


0% 


±9.2% 


-2.9% 


±22% 


-22% 


±35% 




+5% 


±9.6% 


+11% 


±22% 


+30% 


±39% 


0.95 X Mdust 


0% 


±9.6% 


+5.8% 


±22% 


+25% 


±39% 



Note. — Accuracy of the empirical recipes for each model parameter (Ap- 
pendix |XJ and comparison with the results of Draine et al. (2007) f Appendix iBt . 
(1): Each one of the dust model parameters. (2): Offset and rms of the empirical 
fits. By definition, the offset is zero for all parameters, since the corresponding 
estimators are derived from fits applied to our data. The only exception is the 
dust mass, for which the empirical recipe is directly obtained by combining the 
estimators for (U) and Ldust- This leads to dust masses 5% larger than the ones 
yielded by the model (see Appendix IA.4t . If desired, this can be accounted for 
by multiplying by 0.95 (last row). (3): Comparison between the values published 
by D07, obtained by fitting the global SEDs including data from 2MASS, IRAC, 
MIPS and IRAS, and the ones we estimate by applying our recipes to their IRAC 
and MIPS data alone. The first number indicates the relative offset (our value/D07 
value) and the second one the rms. (4): Same as (3), but for the galaxies for which 
submillimcter data from SCUBA were added to the data-set used in D07. 
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Fig. 1. — (a): Radial extinction profiles for the SINGS galaxies, computed from the TIR-to- 
FUV ratio. The horizontal axis is normalized to the optical size of each galaxy. The dashed 
curves show the median profiles in each bin of morphological types, when enough galaxies 
arc available. Top row: Extinction computed Tising the fit of Buat et al. (2005), which 
is independent of the SFH. Bottom row: Extinction computed from the SFH-dependent 
calibration of Cortese et al. (2008). The conversion law is explicitly derived for each annulus 
depending on its observed (FUV — 3.6 fim) color, in order to take into account the dust- 
heating contributed by old stars. For clarity, the profiles of the edge on galaxies NGC 5866 
and NGC 3190 have been offsetted by — Imag in the top panels (but not in the bottom 
ones) . 
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Fig. 2. — Radial gradients of the attenuation in the FUV, measured by applying a linear fit 
to the ApuY profiles, after excluding the bulges. The gradients are expressed in mag/kpc (a), 
in units of the R25 radius (b), and in units of the radial exponential scale-length of the stellar 
mass distribution (c). The latter was derived from the profiles measured on the degraded 
3.6 //m images (bulges excluded). The overall distribution of gradients for all Hubble types 
is shown in the histograms at the right, along with the corresponding median values. For 
comparison, the solid and dotted lines in each panel show the mean and 1-a limits of the 
distribution of attenuation gradients derived by Muiioz-Mateos et al. (2007) for a sample of 
161 nearby, face on spirals. In that paper we relied on the (FUV — NUV) profiles as indirect 
tracers of the extinction, and 2MASS X-band ones as a proxy for the stellar mass. 
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Fig. 3. — Ratio of the total infrared to FUV luminosity as a function of the UV spectral 
slope (/3glx) and the UV color, for all the profiles in the sample. The rightmost y-axis shows 
the internal extinction in the FUV for a given TIR-to-FUV ratio, according to the fits of 
Buat et al. (2005), only valid for late- type systems (see text), (a): Galaxies are sorted out 
into different morphological types. For the disk-like galaxies, a color scheme is used to map 
the intrinsic (FUV — 3.6/im) color, corrected for internal extinction using the age-dependent 
calibration of Cortese et al. (2008). The solid hne corresponds to a fit to the data-points 
bluer than (FUV— NUV) =0.9, excluding three starburst galaxies (see text). The dot-dashed 
line shows the relation found by Boissier et al. (2007) for spiral galaxies using GALEX and 
IRAS data. The dashed line is the mean relation for the starburst galaxies of Meurer et 
al. (1999). The three dotted lines are model predictions by Kong et al. (2004) for different 
values of the birthrate parameter b. The mean error bars (including zero-point errors in 
all bands) are also shown, (b): Individual radial tracks followed by some galaxies in the 
diagram: an early-type spiral (NGC 3031), a late-type one (NGC 2403), an anemic spiral 
(NGC 4826; van den Bergh 1976), and two starbursts (NGC 4536 and NGC 4631, the latter 
being edge-on). The innermost point of each profile is marked for reference. The integrated 
colors of each galaxy measured by Dale et al. (2007) are shown with stars. 
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Fig. 4. — Relation between the intrinsic (FUV — 3.6 /im) color (i.e. corrected for internal 
extinction) and the perpendicular distance to the mean relation for starbursts in the IRX- 
P diagram (see Fig. [3]). Negative distances correspond to data-points below the starbursts 
curve. Different symbols and colors are used depending on the Hubble type and the observed 
UV color. 
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Fig. 5. — (a): Oxygen abundances as a function of the TIR-to-FUV ratio for the radial 
profiles of 22 SINGS galaxies with metallicity gradients in Moustakas et al. (2009, in prepa- 
ration). Data-points lying outside the spatial range covered by the HII regions used to 
compute 0/H (typically beyond the optical size) have not been plotted, as their extrapo- 
lated metallicities are not reliable. Filled circles correspond to moderately face-on galaxies 
{b/a > 0.5), whereas open diamonds indicate those being more edge-on {b/a < 0.5). The 
solid line is a fit to the values of moderately face-on galaxies. The empirical relations of 
Cortese et al. (2006), Boissier et al. (2007) and Heckman et al. (1998) — this latter one valid 
for starbursts alone — are also shown, (b): Relation between metallicity and the observed 
UV color. The symbol code is the same as in panel (a). The solid curve is not a fit to the 
data in this plot; it is the result of the combination of the fit shown in panel (a) with the 
IRX-/? relation in Eq. [21 Median errobars (including zero-point errors) are shown in each 
panel. Fits to the starburst galaxies of Calzetti et al. (1994) and Heckman et al. (1998) are 
also plotted. 
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Fig. 6. — Observed SEDs and model spectra at different galactocentric distances, witfiout 
correcting for inclination. A sample plot for NGC 7331 is shown; additional figures for 
the remaining galaxies can be found in the electronic version of the journal. The observed 
flux densities are marked with rectangular boxes, their width showing the corresponding 
bandpass and their height showing the observational errors. The best-fitting model spectra 
(solid line) results from the combination of a diffuse component (dashed line), the emission 
from hot dust (dash-dotted line) and the stellar emission (dotted line). The filled circles 
correspond to the model spectra convolved with each bandpass. 
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Fig. 7. — Radial distribution of the different parameters derived from the dust models. A 
sample plot is shown for NGC 7331; additional figures for the remaining galaxies can be 
found in the electronic version of the journal. The different fines show the best-fitting values 
at each radius, and the gray bands show the estimated uncertainties, (a): Fraction of the 
dust mass in the form of PAHs. Regions in which the stellar emission at 8 /im is more than 
twice the one from PAHs are marked (when present) with a dotted fine instead of a sofid one. 
The FWHM of the MIPS 160 /xm band (38") is marked with an horizontal segment along 
the major axis, (b): Fraction 7 of the total dust mass heated by very intense starlight. The 
right vertical axis shows the fraction of the dust luminosity contributed by regions exposed 
to radiation fields with U > ?7min and U > 10^ (their errors are not shown for clarity), (c): 
Scale factor of the minimum (t/min) and average {{U)) starlight intensity heating the dust, in 
units of the local MW radiation field. The errors of {U) are not shown, (d): Dust mass (left) 
and luminosity (right) surface densities, both corrected for inclination. The open circles 
show the dust luminosity profiles obtained with the photometric estimator of DL07. 
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Fig. 8. — Radial profiles of the different parameters involved in the dust models. Galaxies 
have been grouped according to their morphological type. The different parameters are the 
PAH abundance (gpAn), the fraction of the dust mass exposed to intense radiation fields 
(7), the intensity of the 'diffuse' heating starlight, in units of the local MW radiation field 

(f^iin)j the deprojected dust mass surface density (J^Mdust) ^^^^ light-to-mass ratio 

(Ldust/Miust)- The radial coordinate is normalized to the optical size of each galaxy. The 
thick dashed lines show the median profiles in each panel. No median profiles are shown for 
ellipticals and lenticulars, due to the small number of galaxies in that bin. The median Umin 
profile for Sdm-Im galaxies is above most of the profiles displayed in that panel due to some 
galaxies having Umin > 5, which are not shown here for clarity. 
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Fig. 9. — Radial exponential scale-lengths of the dust profiles for different morphological 
types. The scale-lengths were derived from linear fits applied to the profiles after having 
excluded the bulges (see text). They are expressed in kpc (a), in units of the optical radius 
(b) and in units of the stellar scale- length (c), measured on the 3.6 //.m profiles derived from 
the degraded IRAC images (bulges excluded). The global distributions for the whole sample, 
as well as the corresponding median values, are shown to the right. 
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Fig. 10. — Attenuation in the FUV as a function of the projected dust mass surface density. 
Black points in each panel show the profiles of galaxies with the corresponding Hubble type. 
The innermost point of each profile is surrounded by a larger circle. As a reference, the 
data for the whole set of annular regions are shown as gray points in each of the panels. 
The dashed line is a fit of the lower envelope of the whole data-set, performed with the 
boundary-fitting code presented in Cardiel et al. (2009) (see text for details). 
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Fig. 11. — Attenuation in the FUV as a function of the /j,(70/im)//j,(160 ;um) color. Black 
points in each panel show the profiles of galaxies with the corresponding Hubble type. The 
innermost point of each profile is surrounded by a larger circle. As a reference, the data for 
the whole set of annular regions are shown as gray points in each of the panels. 
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Fig. 12. — Fraction of the dust mass contributed by PAHs against the metalhcity at each ra- 
dial bin. Data-points with gpAH = 4.6% are 'saturated', meaning that their PAH abundances 
correspond to the current upper hmit of the models (see Section [4. 3. 4p . For these points we 
have recomputed gpAH from the empirical fit presented in Appendix IA.2I (see also Fig. [T8|) . 
The dashed line is a bisector linear fit applied to the data-points below 12 + log{0/H) = 9, 
and corresponds to gpAH oc (0/H)° '^^°''^. Four galaxies for which the gpAH vs. (0/H) trend 
is reversed have been highlighted (see text for details). 
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Fig. 13. — PAH abundance against tlie observed (a) and extinction-corrected (b) (FUV — 
3.6 fj,m) color. Data-points are sorted out depending on their surface brightness in the FUV 
band (in magarcsec"^). To allow for a fair comparison between both panels, an offset of 
1.5 mag has been applied to the surface brightness value chosen to classify the data-points, 
since that is the typical FUV attenuation found in our sample. 
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Fig. 14. — PAH abundance against the photometric estimator of the 8 fim equivalent width, 
computed following the prescriptions of Engelbracht et al. (2008). Different symbols corre- 
spond to different fractions of the total dust mass heated by very intense starlight. 
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Fig. 15. — Radial variation of the dust-to-gas ratio for the SINGS-THINGS galaxies, divided 
into different bins of morphological types. The radial axis is normalized to the optical size 
of each galaxy, as given by the R25 radius. Galaxies with available HI and CO profiles are 
shown with solid lines, whereas those lacking CO data appear as dashed lines (their profiles 
are then upper limits, although most of the molecular gas is expected to be located in the 
innermost regions). In all but two galaxies a metallicity-dcpcndent CO-to-H2 conversion 
factor changing with radius was used (see text). In all plots throughout this paper explicitly 
involving (0/H), we only trust metallicity values within a restricted spatial range defined 
by the HII regions used to determine the metallicity gradients. In this figure, however, the 
accuracy of the oxygen abundance is less critical, since it constitutes a second-order effect in 
comparison with the assumption of a constant CO-to-H2 factor. Therefore, in this particular 
plot we have not applied these restrictions to the metallicity values when computing the H2 
mass surface densities. 
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Fig. 16. — Dust-to-gas ratio for the SINGS-THINGS galaxies against the metallicity at 
different galactocentric distances. Data-points belonging to the same galaxy have been 
connected. Different symbol shapes correspond to different morphological types. The dashed 
line marks the simple scaling law described by Eq. [6l whereas the dotted lines show variations 
of a factor of 2 around it. The dot-dashed line results from a linear fit to the data, (a): 
Galaxies for which both HI and CO data are available, (b): Objects lacking CO data. The 
values at the innermost (i.e. more metallic) regions are upper limits. 



-69- 




0.01 0.1 1 0.01 0.1 1 

P24-0.14P8.0 P24-0.14P8.0 

Fig. 17. — (a): Fraction of the total dust luminosity contributed by grains in regions where 
U > 10^ as a function of the observable P24 — O.MPgo (see text). The solid line is the 
theoretical fitting function provided by Draine & Li (2007). (b): Relation between the 
fraction of the dust mass enclosed in high-intensity regions and P24 ~ O.MPg.o- The solid line 
is the best fit-function of our radial profiles. 
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Fig. 18. — (a): Fraction of the dust mass contributed by PAHs as a function of Pg.o (see 
text). The data-points have been divided into three different bins, according to their value 
of P24 — 0.14P8.0; which is a proxy for the fraction of dust mass heated by intense starhght. 
(b): PAH abundance as a function of a hnear combination of Pg.o and P24, resulting from 
a Principal Component Analysis of the data. The best-fit hne is also shown. Data-points 
whose values of Qpah saturate at the lower and upper limits of the dust model (0.4% and 
4.6%, respectively) were not taken into account when performing the fit. 
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Fig. 19. — Scale factor of the dust-weighted mean starhght intensity —in units of the local 
MW radiation field— against the ratio of the luminosities at 70 and 160 /xm. The best fitting 
hne is also shown. Data-points with {U) ^ 0.7 have been excluded from the fit (see text for 
details) . 
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Fig. 20. — Comparison of the dust mass surface densities obtained from the model- fitting 
technique and from the estimator given in Appendix Dl either with the IRAC+MIPS 
cahbration (Eq. I A 71 panel a) or with the MlPS-only one (Eq. IA81 panel b). Note that 
since we are dealing with surface densities, the results do not depend on the distance to the 
sources. Regions for which the fitting procedure 'saturates' at the lowest value of the diffuse 
starlight intensity ?7min have been marked. 
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Fig. 21. — Comparison of the values of the dust model parameters derived by D07 with the 
ones estimated using our recipes. The former were obtained by using global photometry from 
2MASS, IRAS and Spitzer to constrain the dust models. Our estimates were computed by 
applying the empirical formulae described in Appendix |X] to the integrated IRAC and MIPS 
data alone. Open circles correspond to galaxies for which submillimeter data from SCUBA 
were used by DOT. The mean offsets and scatter are quoted in Table (a): Abundance 
of PAHs. NGC 1291, NGC 1316, NGC 4125 and NGC 4594 have been excluded, since 
their IRAC 8.0 fim flux densities are highly contaminated by stellar emission, leading to very 
uncertain values of gpAH- (b): Fraction of the dust mass heated by intense starlight, (c): 
Dust-weighted mean starlight intensity, in units of the local MW radiation field, (d): Total 
dust mass. 



